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TACTICAL  TARGET  ACQUISITION  MODEL  (TATAC) 


I •  INTRODUCTION 


This  volume  documents  the  operating  procedures,  input  requirements, 
and  output  formats  for  the  Tactical  Target  Acquisition  (TATAC)  computer 
model.  After  defining  the  problem  this  guide  can  be  used  to  assist  the 
user  in  setting  up  the  necessary  steps  to  make  a  computer  run. 

The  computer  program  is  comprised  of  the  following  sensor  system 
model s : 


Visual  Observer  (VISOB) 

Forward-Looking  Infrared  (FLIR) 

Television  (TV) 

-  Active  (Illuminated) 

-  Passive  (Daylight) 

Forward-Looking  Radar  (FLR) 

-  Moving  Target  Indicator  (MTI)  Mode 

-  Non-MTI  Mode 
Synthetic  Aperture  Radar 

Section  II  of  this  guide  describes  the  operation  and  structure  of  the 
program  in  terms  of: 

Execution  List  (Input) 

Output 
Library  Oata 
Fixed  Data 

Section  III  describes  the  logical  structure  of  the  programs.  Section 
IV  contains  a  listing  of  the  fixed  data.  Section  V  contains  a  listing  of 
the  library  data.  Section  VI  contains  sample  problems  including  the  exe¬ 
cution  list  and  output  from  the  computer  run  for  each  problem  to  demon¬ 
strate  the  method  of  using  the  model.  Section  VII  contains  a  listing  of 
all  programs  in  the  model. 
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II.  PROGRAM  OPERATION 


This  section  describes  the  procedures  necessary  to  run  the  TATAC  com¬ 
puter  model.  Detailed  in  this  section  are; 

Description  and  form  of  the  execution  list; 

Form  of  output  generated  by  a  computer  run; 

Library  data  description;  and 
Fixed  data  description. 

Figure  1  is  a  diagram  of  the  program  deck.  Figure  2  summarizes  the 
Information  flow  in  the  program. 


Figure  1.  TATAC  Deck  Structure 
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Figure  2.  TATAC  Program  Information  Flow  (Continued) 


A.  EXECUTION  LIST 

The  execution  list  follows  •she  program  deck  as  shown  in  Figure  1. 

It  begins  with  an  EXECUTE  card  and  ends  with  an  ENDRUN  card.  This  list 
is  the  only  portion  of  the  data  which  must  be  formulated  by  the  user 
before  each  computer  run.  As  a  reference  to  the  following  description, 
a  sample  execution  list  is  shown  in  Figure  3. 

The  execution  list  performs  the  following  functions  in  the  model: 

Determines  the  type  of  run  to  be  executed;  and 
Denotes  all  data  of  the  library  to  be  modified. 

These  functions  are  described  in  the  following  paragraphs. 


The  first  card  in  the  execution  list  is  the  EXECUTE  card  (Figure 
3,  card  1).  Besides  indicating  to  the  computer  model  that  this  is  the 
beginning  of  the  execution  list,  this  card  conveys  the  type  of  run  to  be 
executed.  The  following  sensor  runs  may  be  executed: 

Visual  Observer  (1); 

Forward-Looking  Infrared  (2); 

Active  (illuminated)  Television  (3); 

Passive  (daylight)  Television  (4); 

Forward-Looking  Radar,  MTI  Mode  (5); 

Forward-Looking  Radar,  non-MTI  Mode  (6);  and 

Synthetic  Aperture  Radar  (7). 

To  indicate  the  sensor  run  being  set  up,  the  user  places  one  of 
the  single-integer  codes,  shown  in  parenthesis  in  the  above  list,  into 
column  11  of  the  EXECUTE  card. 

2.  Data  Entry 

Besides  the  EXECUTE  and  ENDRUN  cards  the  only  other  types  of  cards 
that  may  be  needed  are  "data  name"  and  "data  modification"  cards.  A 
"data  mod: fi cation"  card(s)  is  used  in  conjunction  with  a  "data  name"  card. 
That  is,  a  "data  name"  card  need  never  be  used  unless  a  "data  modification" 
card(s)  immediately  follows  it. 

a.  Data  Name  Cards 

The  "data  name"  card  (Figure  3,  Card  2)  has  the  legal  name  of 
the  data  being  modified  in  columns  1-15  (left  justified).  Legal  names 
for  data  used  in  the  program  are  given  in  Table  1.  The  user  then  places 
an  integer  in  columns  20-30  (right  justified).  This  integer  tells  the 
input  program  how  many  modifications  are  to  be  made  to  that  particular 
library  entry.  Then  for  each  modification  a  "data  modification"  card 
follows.  Recall  that  a  "data  name"  card  is  only  used  if  one  or  more 
modifications  are  made  to  that  particular  library  "data  name." 


TABLE  1.  LEGAL  DATA  NAMES  * 


OPERATIONAL  VAR 
TARGET 
BACKGROUND 
ENVIRONMENTAL  ** 
TERRAIN  ** 

SEARCH 

SCNSOR 


Only  the  first  four  characters  are  checked 
by  the  program.  Thus,  all  data  names  in 
the  list  may  be  shortened. 

Environmental  and  terrain  data  are  com¬ 
bined  into  one  array  -  either  name  may  be 
used  as  a  legal  data  name  for  this  array. 


b.  Data  Modification  Cards 


Each  of  these  cards  (Figure  3,  cards  3-4)  contains 
the  number  of  the  variable  to  be  modified  in  columns  9-10  (right 
justified)  and  the  new  value  assigned  to  it  in  columns  11  -  20  in  E10.4 
format  (right  justified).  The  variable  number  and  its  type  are  found 
in  the  library  data  tables  in  Subsection  C  of  this  Section.  The  number 
of  "data  modification"  cards  which  follow  a  "data'name"  card  must  agree 
with  the  integer  placed  in  columns  29  -  30  of  the  "data  name"  card. 

3.  Setting  Up  an  Execution  List 

This  subsection  lists  the  steps  which  the  user  should  follow 
when  setting  up  an  execution  list  to  run  the  TATAC  computer  model. 


Defining  the  problem: 


(1)  Define  the  problem  to  be  run  including  the  sensor  sys¬ 
tem,  operational  variables,  target,  background,  environ 
mental/terrain,  and  search  variables. 


(2) 


Determine  if  the  data  needed  is  in  the  library, 
"data  name"  and  "data  modification"  cards  for  new 
entries  will  be  required. 


If  not 


Filling  out  the  execution  list: 


(3)  List  the  EXECUTE  card  including  the  sensor  number  in 
column  11  for  the  run  to  be  made. 

(4)  Insert  all  "data  name"  cards  along  with  the  proper 
number  of  "data  modification"  cards  for  each. 

(5)  List  the  ENDRUN  card. 


B.  OUTPUT 

There  are  basically  three  categories  of  output  from  the  computer 

model : 


List  of  data  modifications; 

Error  messages;  and 
Standard  model  output. 

Types  of  each  output  applicable  to  a  specific  run  will  be  auto¬ 
matically  generated  by  the  program.  A  discussion  of  the  categories 
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of  output  will  be  given  in  the  next  three  subsections.  Following  these 
discussions  will  be  a  discussion  on  variations  of  standard  model  output 
among  the  sensor  models.  Reference  will  be  made  at  that  point  to  sample 
problems  which  illustrate  the  distinctions  in  output  generated  by  com¬ 
puter  runs. 

1.  List  of  Data  Modifications 

Preceding  all  standard  output,  a  listing  of  data  modifications 
made  through  the  use  of  "data  modification"  cards  in  the  execution  list 
will  be  given.  If  no  data  modifications  have  been  made,  no  output  of 
this  category  will  be  generated.  The  list  will  include  the  legal  data 
name  for  which  modifications  were  made,  the  number  of  the  variable 
modified,  IVAR  (found  in  the  library  data  table  for  the  applicable  data 
name),  and  the  value  assigned  to  that  variable.  Space  for  output  is 
reserved  for  as  many  as  ten  modifications  per  data  name  used.  If  more 
than  ten  "data  modification"  cards  follow  a  "data  name"  card, a  message 
is  printed  alerting  the  user  to  this  occurrence.  This  in  no  way  affects 
the  program  run.  It  simply  means  that  all  changes  made  after  the  tenth 
one,  will  not  be  listed  on  the  output.  There  are  two  ways  to  avoid  this 
situation.  One  is  to  repunch  the  entire  portion  of  the  library  data 
array  applicable  to  the  library  data  name  in  question.  The  other  is  to 
distribute  the  "data  modification"  cards  behind  multiple  "data  name" 
cards  containing  the  same  data  name.  For  example:  if  15  data  modifi¬ 
cations  are  to  be  made  to  a  library  entry,  set  up  a  "data  name"  card 
with  a  "10"  in  columns  29  -  30,  followed  by  ten  "data  modification" 
cards  for  the  first  10  modifications  and  place  the  remaining  five  behind 
another  "data  name"  card  with  a  "5"  in  columns  29  -  30. 

2.  Error  Messages 

At  a  number  of  points  in  the  program,  checks  are  made  on  various 
input  data  to  determine  if  it  is  in  the  correct  form.  If  there  is  an 
error  the  program  will  print  out  a  message.  Generally,  the  program  will 
attempt  to  assign  a  default  value  so  as  to  continue  with  the  run. 

These  messages  may  appear  anywhere  in  the  output  and  usually  give  an 
indication  of  where  the  problem  occurred.  These  messages  are  listed 
below,  followed  by  an  explanation  of  the  cause. 

EXECUTE  CARD  MISSING  OR  OUT  OF  ORDER 

The  first  card  of  the  execution  list  must  be  the  EXECUTE 
card  which  includes  the  sensor  number.  If  the  first  card 
read  is  not  the  EXECUTE  card  then  the  above  message  is 
printed  and  the  sensor  number  defaults  to  the  value  "1". 

INVALID  SENSOR  NUMBER  -  (value) 


Valid  sensor  numbers  range  from  "1"  to  "7".  Any  sensor 
number  input  outside  this  range  will  result  in  the  above 
message.  The  invalid  sensor  number  read  is  printed.  The 
sensor  number  defaults  to  "1". 
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INVALID  OATA  NAME  -  (name) 

The  name  read  does  not  match  one  of  the  legal  data 
names  given  in  Table  1.  The  invalid  name  read  is 
printed. 

VARIABLE  NUMBER  OUTSIDE  RANGE  OF  DATA  NAME  ARRAY  -  (value) 

The  variable  number  in  columns  9  -  10  of  a  "data  modi¬ 
fication"  card  is  outside  the  range  of  the  array 
reserved  for  the  data  name  appearing  on  the  previous 
"data  name"  card.  The  read  value  is  printed. 

* 

MORE  MODIFICATIONS  THAN  LISTED 

This  refers  only  to  the  data  modification  "listing."  It 
does  not  imply  an  error  in  the  run.  Space  is  reserved 
to  print  out  a  maximum  of  10  data  modifications  per  "data 
•  name"  card.  If  more  than  10  "data  modification"  cards 
are  required  for  one  "data  name"  card,  the  user  should 
place  the  remaining  "data  modification"  cards  after 
another  "data  rame"  card  containing  the  same  data  name 
as  for  the  first  10  modifications.  For  many  modifica¬ 
tions  it  may  be  more  efficient  to  repunch  that  section 
of  the  library. 

DEPRESSION  ANGLE  TCO  STEEP,  PHID  RESET  =  (value) 

This  message  is  applicable  to  all  non-radar,  fixed  sensor 
depression  angle  runs.  If  the  input  depression  angle, 
PHID,  is  too  large,  the  computed  dive  profile  point  may 
be  reached  before  the  target  enters  the  sensor  footprint. 
The  program  will  then  set  the  angle  back  to  where  the 
target  just  enters  the  footprint  when  the  dive  profile 
point  is  reached.  The  new  angle,  in  degrees,  is  printed. 

OFFSET  (Y)  GREATER  THAN  YMAX.  SET  Y  =  (value) 

This  message  is  applicable  to  all  models  using  a  fixed 
sensor  depression  angle  with  the  exception  of  the  Syn¬ 
thetic  Aperture  Radar  model.  For  a  given  depression 
angle,  horizontal  beamwidth,  and  vertical  beainwidth,  a 
ground  sensor  footprint  is  calculated.  The  cross-track 
distance  of  the  footprint  at  the  leading  edge  is  calcu¬ 
lated  as  YMAX .  If  the  input  offset  (Y)  is  greater  than 
YMAX,  the  target  never  posses  into  the  ground  footprint. 
If  such  a  condition  is  encountered,  Y  defaults  to  the 
value  printed. 
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ONLY  POSITIVE  OFFSETS  CONSIDERED,  USED  ABS.  VALUE 

Offsets  are  considered  to  be  symmetric  with  respect  to 
flight  path.  Therefore,  if  a  negative  offset  is  input 
the  program  will  convert  it  to  its  absolute  value. 

TUBE  SATURATION,  NEW  FNUM  «  (value) 

Illumination  level  great  enough  to  saturate  tube.  Lens 
is  stopped  down  and  another  trial  is  made  to  see  if  the 
Illumination  is  within  tube  operating  limits. 

THRESHOLD  SPEED  IS  GREATER  THAN  1/2  BLIND  SPEED 

The  MTI  filter  threshold  velocity  exceeds  one-half  the 
calculated  value  of  blind  speed.  This  means  that  the 
filter  has  no  passband  and  that  detection  is  impossible. 

AN  MF  VALUE  WAS  NOT  INPUT  SO  MF  IS  SET  -  12 

For  Forward-Looking  Radar  (MTI  mode)  a  filter  rolloff 
value,  MF,  was  not  input  (or  input  as  zero).  The 
program  defaults  to  a  value  of  12. 

Y  TOO  SMALL,  INCREASED  TO  BE  1.5  SWATH  WIDTHS 

The  offset,  Y,  is  measured  to  the  far  edge  of  the  swath 
width.  No  information  is  obtained  from  the  area  below 
or  close  to  the  aircraft.  For  offsets  less  than  1  -  1/2 
swath  width,  Y  is  reset  to  1.5  times  the  input  swath 
width.  Typical  offset  values  for  SAR  are  on  the  order  of 
10  -  30  nautical  miles. 


3.  Standard  Model  Output 

There  are  basically  three  types  of  standard  model  output  generated 
by  a  computer  run.  All  three  mey  not  be  applicable  to  all  sensor  models. 
The  discussion  below  and  reference  to  Figure  4  should  familiarize  the  user 
with  this  output.  Standard  model  output  follows  the  list  of  data  modifi¬ 
cations  (if  any). 

a.  Variable  Description 

The  first  lines  in  this  output  section  (Figure  4,  lines  12  - 
27)  for  any  run  consist  of  a  brief  description  of  the  variables  which 
head  a  given  column.  Further  description  may  be  found  in  Volume  I 
of  this  report. 
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A*  F«rVAPO-LOO.f  TVf  RADAR.  \JTV-YX  t 
7«  SYVTHFT  IF  AprpxUFF  PAflAp 

<*\*F  I  VAR  VALHF 
"prp  o  .soopf+pa 
opro  4  .poppf*oo 

F  *  *1 V  I  vijv  LA'iMPM  PPTVX 

If  s  D1VF  BFFIVR 

Ills  (PI1FMY  P9S  IT  I  ON  I 

tv  *  XAPFFX  PASSFO  r>nx  OF  F^VCFIXFO  DFPPFSSIDN  ANf-LPI 

V  *  SFAprw  Al.TiTi'nr  aruiM/rn 

VI  *  CLIf'P  TO  ALT I Tl'DF  pep  INS 

X  «  AL^VP-TPACK  PRplINO  PTSTANCF  TO  TARPFT*  f1  T. 

XY  *  FpnijNP  FAvr-F  XO  TARPFX.FT. 

TI^Fr  X  T F-F  orF«RF  LAUNCH, SFC. 

PLPS  =  PRWBAP1LITY  TARfFX  IS  V'TTMIN  I.OS 
PFPV=  PROPAOiLIXY  TARPFT  IS  VITHTN  FOV 
PP  *  SFARpu  XFR"  PP^oadilTTY 
POO  *  0  TSCR  I F I VA®  TLTTYCOFXFPXI  "**i  ) 

P?-R  *  D I SOR  If*  1 N'AP  TL  IXYf  PFO^rv  IT  I" v> 

PAH  s  CM^MLATIVF  Prr°APlLITY  °F  OFTFOTI^V 
PAR  s  PMMJLAXIVF  ppoCAC  )i  txy  af  FFPAPVTXTO.M 


I 

T  I 

III 

TV 

V 

VI 

P7*P. 

OP1  S. 

901  s. 

or  1  5 

.  S070I. 

P7P1 O. 

P7f  P. 

991  S. 

9*1S. 

QP1  5 

.  50701 . 

07010. 

X 

XY 

TIVF 

PI  P  S 

PFO\/ 

po 

pon 

pop 

PAP 

pop 

?Of3. 

P9A0. 

•  OS 

1  .pop 

1  .pop 

1.000  1 

.  OOP 

1 .  ppp 

1  .  POP 

1  . 

nnn 

0770. 

9770. 

n .  aa 

1 .  nnn 

1 .  non 

1 . Pop  1 

.  POP 

.017 

1  .POP 

1  . 

nnn 

900?  • 

o*)91. 

11. 9° 

,  OGE, 

.  907 

.  9P0 

.  1  50 

n.  pon 

.  1  AA 

P. 

nnn 

10AAP. 

1 PAAP. 

IP.  9° 

.  007 

.  000 

.  97  A 

•  PO  A 

0.  POP 

.  1  A° 

P. 

nnn 

1 10A9. 

1 10  A?. 

lO.oo 

.  OOJ 

1  .nor 

.  OA7 

.ns? 

f).  OOP 

.  1  AO 

O. 

noo 

1PP1 A. 

IPP1 A. 

1  A.  oo 

,  OP  o 

l .  non 

.  OSP 

.  00  0 

n.  pon 

.  1  A  9 

0. 

non 

1PA91  . 

IPAO^  . 

1  S.  oc 

.  OP  A 

1 .  non 

.  QA7 

.017 

n.  000 

.  1  AO 

0. 

nnn 

iooas. 

100AS. 

1  A  .  OF 

.  oc  7 

1  .onn 

.  97  A 

.01  1 

p.  oop 

.  1  AO 

0. 

nnn 

1 AOAP. 

»  APAO. 

1  7  .  OP 

.  orp 

1 .  ppp 

.  99  A 

.  00a 

0.  npo 

.  1  A9 

p. 

nnn 

|A7| A. 

I A7 1  A. 

» «.  op 

.  077 

1  .ppp 

.OM 

.  OPA 

P.  POP 

.  1  AO 

0. 

nnn 

Figure  4.  Sample  Computer  Output 
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b.  Flight  Profile  Points 

The  first  two  rows  of  numbers  following. the  variable  descrip¬ 
tion  (Figure  4,  lines  29  -  30)  display  profile  points  along  the  flight 
path  for  all  computer  runs  except  the  Synthetic  Aperture  Radar*.  The 
first  of  these  rows  represent  the  along-track  ground  range,  in  feet,  from 
target  to  aircraft  position  for  the  various  profile  points  listed.  The 
second  row  represents  ground  range,  in  feet,  from  target  to  aircraft 
position  for  the  same  profile  points.  In  case  of  zero  offset  these  two 
rows  are  equivalent. 

c.  Ranges,  Time,  and  Probabilities 

Beneath  the  rows  indicating  flight  profile  (or  immediately 
beneath  the  variable  description  for  the  Synthetic  Aperture  run)  is  a 
row  of  variables  (Figure  4,  line  31)  which  head  10  columns  of  numerical 
data.  These  columns  of  numerical  data  represent  a  combination  of  ranges, 
time,  and  probabilities  which  correspond  to  the  variable  description  of 
whichever  variable  heads  a  given  column.  All  the  data  in  one  particular 
row  is  related  in-so-far-as  it  was  all  generated  by  one  iteration  of 
the  program.  Each  can  be  thought  of  as  being  a  function  of  the  range 
given  in  the  first  column. 

4.  Variations  in  Standard  Model  Output 

A  few  special  points  can  be  made  regarding  the  output  described 
in  the  above  sections. 

a.  Non-Radar  Cases 

For  all  non-radar  cases  the  following  notes  apply: 

The  first  row  of  numerical  output  from  the  "ranges, 
time,  and  probabilities"  output  section  was  computed 
at  a  range  which  corresponded  to  one  "glimpse"  before 
launch  point; 

The  second  of  these  rows  of  output  was  computed  at  a 
range  which  corresponds  to  the  first  "glimpse"  after 
the  dive  profile  point;  and 

The  third  of  these  rows  of  output  was  computed  at  a 
range  which  corresponds  to  the  first  "glimpse"  before 
the  dive  profile  point. 


* 


Output  generated  by  a  Synthetic  Aperture  Radar  run  will  be  given 
special  consideration  in  a  following  subsection. 


b.  Forward-Looking  Radar 

As  discussed  in  Volume  I  to  this  report,  all  radar  submodels 
are -essentially  taken  from  f'ARSAM*  with  technology  updates.  Users  familiar 
with  MARSAM  will  recall  that  it  utilized  a  fixed  depression  angle  for  radar 
sensors.  Thus,  for  given  horizontal  and  vertical  beamwidths  a  sensor  foot¬ 
print  on  the  ground  was  clearly  defined.  Results  at  one  range  between  the 
leading  and  trailing  edges  of  this  footprint  were  output.  In  the  TATAC 
version  a  footprint  is  clearly  defined,  however,  output  will  be  generated 
for. Intermittent  ranges  from  the  trailing  edge  up  to  the  leading  edge. 

The  flight  profile  is  straight  and  level  throughout  and  no  dive  is 
assumed.  Profile  points  normally  printed  thus  appear  as  blanks  in  the 
dive  phase.  Refer  to  Problem  B,  Section  VI  for  an  example  of  output  from 
the  Forward-Looking  Radar  models. 

c.  Synthetic  Aperture  Radar 

The  Synthetic  Aperture  Radar  is  also  a  version  of  a  MARSAM 
submodel.  Again,  there  is  no  dive  profile.  Furthermore,  the  Synthetic 
Aperture  Radar  processes  the  doppler  phase  history  of  the  illuminated 
ground  area  over  the  entire  illumination  period  to  generate  the  display. 

All  returns  from  the  target  are  integrated  into  this  display.  Thus, 
only  one  line  of  output  is  necessary  to  present  the  results.  For  an 
example,  see  Problem  C,  Section  VI. 


C.  PROGRAM  LIBRARY  DATA 

The  program  library  consists  of  a  group  of  DATA  declaration  statements 
found  in  subroutine  INPUT!.  It  contains  information  needed  by  the  program 
to  describe  sensors,  target,  background,  etc.  The  library  may  contain 
more  information  than  is  needed  by  the  program  for  a  particular  run.  A 
complete  list  of  the  library  data  supplied  with  the  program  is  given  in 
Section  V.  Tables  2  through  13  following  this  discussion  define  the 
sets  of  entries  allowed  in  the  library.  These  sets  of  library  data  are: 

Operational  variables; 

Target  data; 

Background  data; 

Environmental /terra in  data; 

Search  data; 

Visual  Observer  sensor  data; 

Forward-Looking  Infrared  sensor  data; 

Television  sensor  data  (active); 


*  "Multiple  Airborne  Reconnaissance  Sensor  Assessment  Model",  Honeywell 
Technical  Report  ASD-TR-G8-3,  February  1968,  Unclassified. 
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Television  sensor  data  (passive); 

Forward-Looking  Radar  (MTI)  sensor  data; 

Forward-Looking  Radar  (non-MTI)  sensor  data;  and 
Synthetic  Aperture  Radar  sensor  data. 

Each  of  these  tables  lists  the  name  of  the  type  data  which  it  des¬ 
cribes,  the  number  of  integers  in  the  data  (MI),  and  the  number  of 
floating-point  numbers  (NF).  For  each  integer  and  floating-point  number 
used  to  define  the  entry,  the  table  lists  the  type  of  data,  the  variable 
number,  the  FORTRAN  symbol  which  is  to  represent  the  variable,  the  units 
in  which  the  data  is  expressed,  and  a  description  of  the  entry. 

TABLE  2.  OPERATIONAL  VARIABLES 
NAME:  OPERATIONAL  VAR:  NI  -  1;  NF  -  14 


VARIABLE 

NUM8ER 


VARIA3LE 

TYPE 

UNITS 

ID 

FD 

T 

FEET 

FD 

HL 

FEET 

FD 

HP 

FEET 

FD 

SFD 

KNOTS 

FD 

SFL 

KNOTS 

FD 

SPC 

KNOTS 

FD 

SPF 

KNOTS 

FU 

SRX 

FEET 

FD 

SR2 

FEET 

FD 

ANCLD 

DEG 

FD 

ANCLC 

DEC 

FD 

FD 

- 

- 

FD 

- 

- 

DESCRIPTION 


UNUSED 

FLIGHT  PATH  OFFSET 
SEARCH  ALTITUDE 
PENETRATION  ALTITUDE 
SPEED  IN  DIVE' 

SPEED  IN  LEVEL  FLIGHT 

SPEED  IN  CLIMB 

SPEED  IN  PENETRATION 

SLANT  RANGE,  MINIMUM  LAUNCH 

SLANT  RANGE,  POP-UP  ATTAINED 

DIVE  ANCLE 

CLIMB  ANGLE 


UNUSED 

UNUSED 


i 


TABLE  3.  TARGET  DATA 
NAME:  TARGET;  NI  -  3;  NF  -  17 


VARIABLE 

NUMBER 

VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

UNITS 

DESCRIPTION 

1 

10 

NJ 

- 

TOTAL  NUMBER  OF  TARGET  ELEMENTS 

2 

ID 

- 

- 

UNUSED 

)  * 

TO 

- 

- 

UNUSED 

4 

TO 

DTE 

FEET 

TARGET  ALONC-TRACK  DIMENSION 

5 

TO 

DTY 

PEET 

TARGET  CROSS-TRACK  DIMENSION 

6 

to 

DTZ 

FEET 

TARGET  HEIGHT 

7 

TO 

RT 

- 

TARGET  REFLECTANCE 

8 

TO 

DELTT 

•k 

TARGET  TEMPERATURE 

9 

TO 

EMMT 

- 

TARGET  EMISSIVITY 

10 

TO 

SIGTX 

MEDIAN  TARGET  RCS  AT  1  GHZ 

U 

TO 

- 

- 

UNUSED 

12 

TO 

- 

UNUSED 

13 

TO 

TSPACE 

FEET 

SPACING  BETWEEN  TARGETS 

14 

TO 

PSI 

DEG 

ANGULAR  DIRECTION  OF  TARGET  RELATIVE  TO  FLIGHT  FATR 

13 

TO 

VTT 

KNOTS 

TARGET  VELOCITY 

16 

TO 

- 

- 

UNUSED 

17 

TO 

- 

UNUSED 

18 

TO 

- 

- 

UNUSED 

19 

TO 

- 

- 

UNUSED 

20 

TO 

UNUSED 

TABLE  4.  BACKGROUND  DATA 
NAME:  BACKGROUND;  NI  -  1 ;  NF  -  6 


VARIABLE 

NUMBER 

VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

— 

UNITS 

DESCRIPTION 

1 

ID 

UNUSED 

2 

TO 

RB 

- 

REFLECTANCE  OF  BACKGROUND 

3 

TO 

DELTB 

*K 

TEMPERATURE  OF  BACKGROUND 

4 

TO 

EKMB 

- 

©IISSIVITY  OF  BACKGROUND 

5 

TO 

- 

- 

UNUSED 

•6 

TO 

WP 

METERS 

EXTENT  OF  PRIMARY  BACKGROUND  AROUND  TARGET 

7 

TO 

- 

UNUSED 
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TABLE  5.  ENVIRONMENTAL/TERRAIN  DATA 
NAME:  ENVIRONMENTAL  OR  TERRAIN;  NI  -  5;  NF  -  8 


VARIABLE 

NUMBER 

VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

UNITS 

DESCRIPTION 

1 

TD 

I CLOUD 

- 

CLOUD  COVERACE  (1  -  CLEAR,  2  -  PARTIAL  OVERCAST 

3  -  SOLID  OVERCAST)  ‘  ' 

2 

ID 

ISUNAN 

- 

SUN  ANCLE  ABOVE  HORIZON  (1  -  90*.  2  -  30* ,  3  -  11.5*) 

3 

ID 

XCFLOS 

SWITCH  FOR  CLOUD  FREE  LOS  (0  -  DO  NOT  CONSIDER,  1  - 
use  tables)  * 

4 

ID 

ITAT 

SWITCH  FOR  ATMOSPHERIC  TRANSMITTANCE  (0  -  USE 

ANALYTIC  APPROXIMATION,  1  -  USE  EMPIRICAL  DATA) 

5 

ID 

XAZQf 

- 

SUN  AZIMUTH,  USED  IF  ITAT  -  1  (1  -  90*,  2  -  0*) 

6 

n 

RTE 

- 

REFLECTANCE  OF  TERRAIN 

7 

FD 

AKCUI 

DEC. 

AVERAGE  MASKING  ANCLE 

• 

FD 

RATIO 

- 

CULTURAL  KASKINC  RATIO 

9 

FD 

VC 

NM 

METEOROLOGICAL  VISIBILITY 

10 

n 

- 

RELATIVE  HUMIDITY 

11 

TD 

c 

- 

MEASURE  OF  SCENE  COMPLEXITY 

12 

FD 

vw 

K/SEC 

APPARENT  RMS  VELOCITY  DUE  TO  WIND 

11 

FD 

- 

• 

TABLE  6.  SEARCH  DATA 
NAME:  SEARCH;  NI  -  1;  NF  -  9 


VARIABLE 

NUMBER 

VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

UNITS 

DESCRIPTION 

1 

ID 

■ 

UNUSED 

2 

FD 

- 

- 

UNUSED 

3 

FD 

DSX 

FEET 

SEARCH  LENGTH* 

4 

FD 

DSY 

FEET 

SEARCH  WIDTH* 

S 

FD 

WLOC 

FEET 

LOC  WIDTH 

6 

FD 

SIGN 

STANDARD  DEVIATION  OF  ALONO-TRACK  ERROR 

7 

FD 

SIGY 

STANDARD  DEVIATION  OF  CROSS-TRACK  ERROR 

8 

FD 

- 

- 

UNUSED 

9 

FD 

- 

- 

UNUSED 

10 

FD 

- 

- 

UNUSED 

*  F.NTER  -1  for  those  variables  if  FOV  is  to  bo  searched.  The  visual  case 
uses  a  unique  geometry  in  search  along  a  LOC  as  described  in  App  A. 
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TABLE  7.  -VISUAL  SENSOR  DATA 
NAME:  VISOB;  NI  -  2;  NF  -  5 


VARIA8LE 

NUMBER 

VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

UNITS 

DESCRIPTION 

1 

ID 

ICMASK 

. 

SWITCH  FOR  COCKPIT  MASKING  (0  -  NO,  1  -  YES)' 

2 

ID 

- 

- 

UNUSED 

3 

FD 

VISMAX 

FEET 

MAXIMUM  VISUAL  SEARCH  RANGE 

A 

FD 

- 

- 

UNUSED 

5 

FD 

- 

- 

UNUSED 

6 

FD 

XLAMDA 

um 

RESPONSIVE  WAVE  LENGTH 

7 

FD 

- 

- 

UNUSED 

TABLE  8.  FORWARD-LOOKING  INFRARED 
NAME:  FUR;  NI  -  2;  NF  -  IQ 


VARIABLE 

NUMBER 

VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

UNITS 

DESCRIPTION 

1 

ID 

ISTYPE 

SWITCH  FOR  DEPRESSION  ANGLE  (1  -  VARIABLE,  Z  -  FIXED) 

2 

ID 

- 

UNUSED 

3 

FD 

ASPECT 

- 

RATIO  OF  THE  VERTICAL  FOV  TO  HORIZONTAL  FOV 

A 

FD 

THETAV 

DEG 

VERTICAL  FOV  (NARROW) 

5 

FD  . 

THETAV 

DEG 

VERTICAL  FOV  (WIDE) 

6 

FD 

PHID 

DEG 

SENSOR  DEPRESSION  ANCLE 

7 

FD 

ZK1 

- 

LINEAR  CONSTANT  TO  KRT  CURVE 

8 

FD 

ZK2 

- 

EXPONENTIAL  CONSTANT  TO  MRT  CURVE 

9 

FD 

“ 

UNUSED 

10 

FD 

UNUSED 

11 

FD 

- 

- 

UNUSED 

12 

FD 

- 

- 

UNUSED 
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TABLE  9.  ACTIVE  TV  SENSOR  DATA 
NAME:  TELEVISION;  NI  -  2;  NF  -  24 


VARIABLE 

DUMBER 

VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

UNITS 

DESCRIPTION 

1 

to 

ISTYPE 

- 

SWITCH  FOR  DEPRESSION  ANCLE  (1- VARIABLE,  2-  FIXED) 

2 

ID 

- 

- 

UNUSED 

3 

FD 

ASPECT 

- 

RATIO  OF  THE  VERTICAL  FOV  TO  HORIZONTAL  FOV 

4 

ID 

THETAV 

DEG 

VERTICAL  FOV  (NARROW) 

5 

FD 

THETAV 

DEG 

VERTICAL  FOV  (WIDE) 

6 

FD 

PHID 

DEC 

SENSOR  DEPRESSION  ANGLE 

7 

FD 

BAND 

HZ 

BANDWIDTH 

8 

FD 

DIAG 

MM 

DIAGONAL  OF  EFFECTIVE  PHOTOSURFACE  AREA 

9 

FD 

mm 

- 

F  -  NUMBER  (NARROW  FOV) 

10 

FD 

FNUM 

- 

F  -  NUMBER  (WIDE  FOV) 

11 

FD 

CAMMAT 

- 

SLOPE  OF  SICNAL  VS/  XRRADIANCE  CHARACTERISTIC 

12 

FD 

TRANSM 

- 

LENS  SYSTEM  TRANSMITTANCE 

13 

FD 

XFC 

FC 

SENSITIVITY  DATA  ILLUMINANCE 

14 

FD 

XI 

AMPS 

SICNAL  AT  Ht 

15 

FD 

X1MAX 

AMPS 

MAXIMUM  SIGNAL  CAPABILITY 

16 

FD 

XIP 

AMPS 

PRE-AMP  NOISE 

17 

FD 

XLAMDA 

RESPONSIVE  WAVE  LENGTH 

18 

FD 

XNL 

- 

5Z  HORIZONTAL  FREQUENCY  RESPONSE 

19 

FD 

XNR 

- 

RASTER  COUNT 

20 

FD 

POUT- 

WATTS 

EFFECTIVE  POWER  OP  THE  ILLUMINATOR 

21 

FD 

CAIN 

- 

BACKSCATTER  CAIN  RELATIVE  TO  ISOTROPIC  ILLUMINATOR 

22 

FD 

CT 

- 

TUBE  CAIN  WN 

23 

FD 

- 

UNUSED 

24 

FD 

- 

- 

UNUSED 

25 

FD 

- 

- 

UNUSED 

26 

FD 

- 

- 

UNUSED 

TABLE  10.  PASSIVE  TV 
NAME:  TELEVISION;  HI  -  2;  NF  -  24 


232291 

wfl 

VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

UNITS 

ID 

ISTYPE 

- 

2 

ID 

- 

- 

3 

FD 

ASPECT 

- 

4 

FD 

THETAV 

DEC 

5 

FD 

THETAV 

DEC 

6 

FD 

PHID 

DEC 

7 

FD 

BAND 

HZ 

8 

FD 

DIAG 

MM 

9 

FD 

FNUM 

- 

10 

FD 

FNUM 

- 

11 

FD 

CAMMAT 

- 

12 

ID 

TRANSM 

j- 

13 

FD 

1  ! 

XFC 

(PC 

14 

FD  1 

xi  ! 

AMPS 

15  j 

FD 

XIMAX 

AMPS 

16 

FD 

XIP 

AMPS 

17 

FD 

XLAMDA 

ura 

18 

FD 

XML 

- 

19 

FD 

XSR 

- 

20 

FD 

POUT 

WATTS 

21 

FD 

CAIN 

- 

22 

FD 

GT 

- 

23 

FD 

- 

w 

24 

FD 

- 

25 

FD 

- 

- 

26 

FD 

m 

- 

DESCRIPTION 


SWITCH  FOR  DEPRESSION  ANGLE  (1-VARIABLE,  2-  FIXED)) 
UNUSED 

RATIO  OF  THE  VERTICAL  FOV  TO  HORIZONTAL  FOV 
VERTICAL  FOV  (NARROW) 

VERTICAL  FOV  (WIDE) 

SENSOR  DEPRESSION  ANCLE 
BANDWIDTH 

DIAGONAL  OF  EFFECTIVE  PHOTOSURFACE  AREA 
F  -  NUMBER  (NARROW  FOV) 

F  -  NUMBER  (WIDE  FOV) 

SLOPE  OF  SICKAL  VS/  IRRADIAKCE  CHARACTERISTIC 
LENS  SYSTEM  TRANSMITTANCE 
SENSITIVITY  DATA  ILLUMINANCE 
SIGNAL  AT  Hi 

MAXIMUM  SIGNAL  CAPABILITY 
PRE-AMP  NOISE 
RESPONSIVE  WAVE  LENCTH 
52  HORIZONTAL  FREQUENCY  RESPONSE 
RASTER  count 

EFFECTIVE  POWER  OF  THE  ILLUMINATOR 
BACKSCATTER,  CAIN  RELATIVE  TO  ISOTROPIC  ILLUMINATOR  WAV 

TJBE  GAIN 
UNUSED 
UNUSED. 

UNUSED 

UNUSED  _ _ _ _ 


LENCTH 


TABLE  11.  FORWARD-LOOKING  RADAR,  MTI 
NAME:  FLR;  NI  -  2;  NF  -  23 


VARIABLE 

NUMBER 


VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

ID 

ISEAS 

ID 

- 

FD 

FXF 

FD 

GO 

FD 

FNR 

FD 

PRF 

FD 

PX 

FD 

FLAMDA 

FD 

WS 

FD 

PHIM 

FD 

THETAH 

FD 

THETAS 

FD 

THETAV 

FD 

TAUPW 

FD 

TAUP1 

FD 

TAUF2 

FD 

VF 

FD 

VPRF 

FD 

- 

FD 

- 

FD 

CAS 

FD 

- 

FD 

TSM 

FD 

- 

FD 

VTH 

db/oct 


DESCRIPTION 


SEASON  (1  -  WINTER,  2  -  SPRING,  3  -  SUMMER,  4  -  FALL) 

UNUSED 

FILTER  ROLLOFF  VALUE 
MAXIMUM  ANTENNA  GAIN 
RECEIVER  NOISE  POWER 
PULSE  REPETITION  RATE 

PEAK  TRANSMITTER  POWER  ' 

WAVELENGTH 

ANTENNA  SCAN  RATE 

ANTENNA  DEPRESSION  ANGLE  TO  TOP  OF  BEAM 
HORIZONTAL  ANTENNA  BEAMWIDTH 

TOTAL  HORIZONTAL  SCAN  ANGLE,  MEASURE  IN  HORIZONTAL 

VERTICAL  ANTENNA  BEAMWIDTH  | 

PULS3WIDTH  ! 

TRANSMITTER  PLUMBING  ATTENUATION  FACTOR 

RECEIVER  PLUMBING  ATTENUATION  FACTOR 

DOPPLER  SPREAD  DUE  TO  CARRIER  FREQUENCY  DEVIATION 

DOPPLER  SPREAD  DUE  TO  PKF  DEVIATION 

UNUSED 

UNUSED 

FRACTION  OF  TIME  ATTENUATION  NOT  EXCEEDED 
UNUSED 

SIGNAL  TO  NOISE  THRESHOLD 
UNUSED 

THRESHOLD  VELOCITY  (RADIAL) 


TABLE  12.  FORWARD-LOOKING  RADAR,  NON-MTI 
NAME:  FLR;  NI  -  2;  NF  -  23 


VARIABLE 

NUMBER 

VARIABLE 

TYPE 

FORTRAN 

SYMBOL 

UNITS 

DESCRIPTION 

1 

n 

I  SEAS 

- 

SEASON  (1  -  WINTER,  2  -  SPRING,  3  -  SUMMER,  4  -  FALL) 

2 

ID 

- 

- 

UNUSED 

3 

TO 

FKF 

4b/occ 

FILTER  ROLLOFF  VALUE 

4 

FD 

CO 

- 

MAXIMUM  ANTENNA  GAIN 

5 

re 

PNR 

WATTS 

RECEIVER  NOISE  POWER 

6 

re 

PRF 

PPS 

PULSE  REPETITION  RATE 

7 

re 

PX 

WATTS 

PEAK  TRANSMITTER  POWER 

8 

re 

FLAMDA 

METERS 

WAVELENGTH 

9 

re 

VS 

DEG/SEC 

ANTENNA  SCAN  RATE 

10 

re 

PHIM 

DEC 

ANTENNA  DEPRESSION  ANCLE  TO  TOP  OF  BEAM 

11 

re 

THETAH 

DEG 

HORIZONTAL  ANTENNA  BEAMWIDTH 

12 

re 

THETAS 

DEG 

TOTAL  HORIZONTAL  SCAN  ANCLE,  MEASURE  IN  HORIZONTAL 

13 

re 

THETAV 

DEC 

PLANE 

VERTICAL  ANTENNA  BEAMWIDTH 

14 

re 

TAUPW 

SEC 

PULSEWIDTH 

15 

re 

TAUP1 

- 

TRANSMITTER  PLUMBING  ATTENUATION  FACTOR 

16 

re 

TAUP2 

- 

RECEIVER  PLUMBING  ATTENUATION  FACTOR 

17 

re 

VP 

CPS 

DOPPLER  SPREAD  DUE  TO  CARRIER  FREQUENCY  DEVIATION 

18 

re 

VPR? 

CPS 

DOPPLER  SPREAD  DUE  TO  PRF  DEVIATION 

19 

re 

- 

- 

UNUSED 

20 

re 

- 

- 

UNUSED 

21 

re 

CAS 

- 

FRACTION  OF  TIME  ATTENUATION  NOT  EXCEEDED 

22 

re 

- 

- 

UNUSED 

23 

re 

TSM 

db 

SIGNAL  TO  NOISE  THRESHOLD 

24 

re 

- 

- 

UNUSED 

23 

re 

VTH 

M/SEC 

THRESHOLD  velocity  (RADIAL) 
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TABLE  13.  SYNTHETIC  APERTURE  RADAR 
NAME:  SAR;  MI  -  2;  NF  -  23 


VARIABLE 

NUMBER 

VARIABLE 

TYPE 

FORTRAN 

SYKSOL 

UNITS 

DESCRIPTION 

1 

ID 

ISEAS 

- 

SEASON  (1-WINTER.  2 -SPRING,  3-SUMMER,  4-FALL) 

2 

ID 

- 

«• 

UNUSED 

3 

FD 

DX 

METERS 

AZIMUTH  GROUND  RESOLUTION 

4 

ID 

CZER 

- 

UNUSED 

S 

FD 

?m 

WATTS 

RECEIVER  NOISE  POWER 

6 

FD 

FRF 

PPS 

PULSE  REPETITION  RATE 

7 

FD 

PX 

WATTS 

PEAK  TRANSMITTER  POWER 

8 

FD 

RC 

- 

CHIRP  RATIO 

9 

FD 

SOT 

db 

SIGNAL  TO  NOISE  THRESHOLD 

10 

FD 

- 

- 

UNUSED 

11 

FD 

FLAMDA 

METERS 

WAVELENGTH 

12 

FD 

- 

- 

UNUSED 

13 

FD 

- 

- 

UNUSED 

14 

FD 

THETAV 

DEG 

ANTENNA  VERTICAL  FOV 

15 

ID 

TAUPW 

SEC 

PULSEWIDTH 

16 

FD 

TAUF1 

- 

TRANSMITTER  PLUMBING  ATTENUATION  FACTOR 

17 

FD 

TAUP2 

- 

RECEIVER  PLUMBING  ATTENUATION  FACTOR 

16 

FD 

- 

- 

UNUSED 

19 

FD 

ASPECT 

- 

20 

FD 

SQ 

DEC 

SQUINT  ANCLE 

21 

FD 

CAS 

- 

FRACTION  OF  TIME  ATTENUATION  MOT  EXCEEDED 

22 

FD 

su 

- 

FIRST  SIDELOBE  LEVEL  WITH  RESPECT  TO  HAINLOBE 

23 

FD 

SL2 

- 

SECOND  SIDELOBE  LEVEL  WITH  RESPECT  TO  HAINLOBE 

24 

FD 

SL3 

- 

THIRD  SIDELOBE  LEVEL  WITH  RESPECT  TO  HAINLOBE 

25 

FD 

VS 

KM 

CROUND  SWATH  WHICH  CAN  BE  PROCESSED  IN  REAL  TIME 

0.  FIXED  DATA 


The  fixed  data  consists  of  a  group  of  DATA  declaration  statements 
found  in  subroutine  IMPUT1.  It  contains  data  which  is  constant  for  all 
senior  models  in  the  program.  It  should  not  be  modified  by  the  user 
unless  the  program  itself  is  modified  or  the  user  wishes  to  change  some 
of  the  contained  data  which  is  discussed  below.  The  only  method  of 
modifying  the  fixed  data  is  to  repunch  the  cards  containing  the  data. 

As  a  reference  for  the  discussion  of  this  subsection,  see  the  Fixed 
Data  List  In  Section  IV. 

The  first  portion  of  the  fixed  data  contains  50  Integers  which  are 
described  in  Table  14.  This  portion  of  the  fixed  data  is  addressed  as 
array  IC  in  the  program. 


TABLE  14.  INTEGERS  IN  FIXED  DATA  ARRAY  IC 


jlNTEGEK  NUMBER  DESCRIPTION  TY?E  OF  DATA*  1 

i 

MEMORY  ALLOCATED  FOR  OPERATIONAL  VARIABLES 

ID 

2 

MEMORY  ALLOCATED  FOR  OPERATIONAL  VARIABLES 

FD 

3 

MEMORY  ALLOCATED  FOR  A  TARGET  ENTRY 

ID 

4 

MEMORY  ALLOCATED  FOR  A  TARGET  ENTRY 

FD 

3 

MEMORY  ALLOCATED  FOR  A  BACKGROUND  ENTRY 

ID 

6 

MEMORY  ALLOCATED  FOR  A  BACKGROUND  ENTRY 

FD 

7 

MEMORY  ALLOCATED  FOR  AN  ENVIRONMENTAL/TERRAIN  ENTRY 

ID 

8  ' 

MEMORY  ALLOCATED  FOR  AN  ENVIRONMENTAL/TERRAIN  ENTRY 

FD 

9 

MEMORY  ALLOCATED  FOR  SEARCH  VARIABLES 

ID 

10 

MEMORY  ALLOCATED  FOR  SEARCH  VARIABLES 

FD 

11 

MEMORY  ALLOCATED  FOR  A  VISUAL  SENSOR  DrTRT 

ID 

12 

MEMORY  ALLOCATED  FOR  A  VISUAL  SENSOR  ENTRY 

FD 

13 

MEMORY  ALLOCATED  FOR  A  FLIR  SENSOR  ENTRY 

ID 

14 

MEMORY  ALLOCATED  FOR  A  FLIR  SENSOR  ENTRY 

FD 

13 

MEMORY  ALLOCATED  FOR  AN  ACTIVE  TV  SENSOR  ENTRY 

ID 

16 

MEMORY  ALLOCATED  FOR  AN  ACTIVE  TV  SENSOR  ENTRY 

FD 

17 

MEMORY  ALLOCATED  FOR  A  PASSIVE  TV  SENSOR  ENTRY 

ID 

18 

MEMORY  ALLOCATED  FOR  A  PASSIVE  TV  SENSOR  ENTRY 

FD 

19 

MEMORY  ALLOCATED  FOR  A  MTI  FLR  SENSOR  ENTRY 

ID 

20 

MEMORY  ALLOCATED  FOR  A  MTI  FLR  SENSOR  ENTRY 

FD 

21 

KBtORY  ALLOCATED  FOR  A  NON-MTI  FLR  SENSOR  ENTRY 

ID 

22 

MD10RY  ALLOCATED  FOR  A  NON-MTI  FLR  SENSOR  ENTRY 

FD 

23 

MB10RY  ALLOCATED  FOR  A  SAR  SENSOR  ENTRY 

ID 

26 

MEMORY  AILOCATED  FOR  A  SAR  SENSOR  ENTRY 

FD 

25-50 

UNUSED  LOCATIONS 

• 

•  io 

INTEGER;  FD  -  FLOAT INC- POINT 

The  second  portion  of  the  fixed  data  contains  50  additional  integers 
which  are  described  in  Table  15.  This  portion  of  the  fixed  data  Is 
addressed  as  array  ICT  in  the  program. 


TABLE  15.  INTEGERS  IN  FIXED  DATA  ARRAY  ICT 
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TABLE  15.  INTEGERS  IN  FIXED  DATA  ARRAY  ICT  (Continued) 


(  INTECER  NUMBER  DESCRIPTION  *  j 

31 

MEMORY  ALLOCATED  FOR  TAT,I,F.  LUMBER  11  OK  FDTAB  ARRAY 

32 

NUMBER  OF  FIRST  INDEPENDENT  VARIABLES  IN  TABLE  11 

33 

NUMBER  OF  SECOND  INDEPENDENT  VARIABLES  IN  TABLE  11 

34 

MQ»RY  ALLOCATED  FOR  TABLE  LUMBER  12  OF  FDTAB  ARRAY 

33 

NUMBER  OF  FIRST  INDEPENDENT  VARIABLES  IN  TABLE  12 

36 

NUMBER  OF  SECOND  INDEPENDENT  VARIABLES  IK  TABLE  12 

37 

MEMORY  ALLOCATED  FOR  TABLE  NUMBER  13  OF  FDTAB  ARRAY 

38 

NUMBER  OF  FIRST  INDEPENDENT  VARIABLES  IN  TABLE  13 

39 

NUMBER  OF  SECOND  INDEPENDENT  VARIABLES  IN  TABLE  13 

40 

MEMORY  ALLOCATED  FOR  TABLE  NUMBER  14  OF  FDTAB  ARRAY 

41 

HUMBER  OF  FIRST  INDEPENDENT  VARIABLES  IN  TABLE  14 

42  ■ 

NUMBER  OF  SECOND  INDEPENDENT  VARIABLES  IN  TABLE  14 

43-30 

UNUSED  LOCATIONS 

*  A  LIST  OF  THE  TABLE  NUMBERS  AND  DESCRIPTION  FOR  THE  FDTAB 

ARRAY  CAM  BE  FOUND  IN  TABLE  16. 

.  The  third  portion  of  the  fixed  data  contains  all  the  tables  used 
In  the  program.  The  original  14  tables  are  described  in  Table  16. 

This  portion  of  the  fixed  data  is  addressed  by  the  array  FDTAB. 

The  floating-point,  third  portion  of  the  fixed  data  contains  all 
the  tables  used  in  the  program.  The  original  14  tables  are  described 
In  Table  16.  This  portion  of  the  fixed  data  is  addressed  by  the  array 
FDTAB. 
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TABLE  16.  TABLES  IN  FIXED  DATA  ARRAY  FDTAB 


TABLE  t.liMtim  TABUS  SIZE  DESCRIPTION 


1 

88 

TABULAR  VALUES  FOR  Ej  (2)  vh«ta 

. 

RjCz)  -j/"  dt 

2 

8 

PROBABILITY  OF  CLOUB-FREE  LOS  VERSUS  ZENITH 

ANCLE  (DEC.)  FOR  ALTITUDES  LESS  TRAN  OR  EQUAL  TO  5000  FT. 

3 

8 

FROSABILITY  OF  CLOUD-FREE  LOS  VERSUS  ZENITH 

ANCLE  (DEC.)  FOR  ALTITUDES  GREATER  THAN  5000  FEET 

4 

IS 

BACKGROUND  RCS  (M~2)  VERSUS  DEPRESSION  ANGLE  (DEC.) 

AND  RADAR  BAND. 

3 

31 

ATTENUATION  RATE  VERSUS  RAIN  RATE  QM/HR)  AID 

RADAR  RAND. 

« 

14 

ALTITUDE  OF  CLOUD  LAYER  (CEIL INC  OR  FLOOR)  VERSUS 

SEASON 

7 

49 

LIQUID  HATER  CONTENT  (g/n3)  VERSUS  SEASON  AMD 

OCTALS  OF  CLOUD  COVER 

8 

49 

PERCENTAGE  FREQUENCY  OF  CLOUD  COVER  VERSUS  SEASON 

AND  OCTALS  OF  CLOUD  COVER. 

9 

39 

FREQUENCY  OF  RAINFALL  RATE  VERSUS  SEASON  AND 

RAINFALL  RATE  QM/HR.)  SLANT  RANCE  FROM  BOTTOM 

OF  CLOUD  COVER  TO  GROUND  GREATER  THAN  OR  EQUAL 

TO  10  KM. 

10 

39 

FREQUENCY  OF  RAINFALL  RATE  VERSUS  SEASON  AND 

RAINFALL  RATE  (MM/HR)  SLANT  RANCE  FROM  EOTTOM 

OF  CLOUD  COVER  TO  GROUND  LESS  THAN  10  KM. 

11 

181 

DIRECTIONAL  PATH  REFLECTANCE  VERSUS  ALTITUDE  (FT.) 

AND  SLANT  RANCE  (FT.)  FOR  90*  AZIMUTH. 

12 

181 

DIRECTIONAL  PATH  REFLECTANCE  VERSUS  ALTITUDE  (Ft) 

AND  SLANT  RANCE  (FT)  FOR  0*  AZIMUTH. 

13 

8 

REFLECTANCE  FACTOR  VERSUS  ZENITH  ANGLE  (DEC.) 

FOR  90*  AZIKDTH. 

14 

8 

REFLECTANCE  FACTOR  VERSUS  ZENITH  ANCLE  (DEC.) 

FOR  0*  AZD1UTU. 
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III.  LOGICAL  STRUCTURE 


This  program  is  arranged  modularly  with  respect  to  the  sensor 
models.  The  control  section  is  independent,  allowing  sensor  models 
to  be  added  or  changed  without  changing  the  whole  program. 

Four  main  data  arrays  form  the  basic  structure  of  the  TATAC  model 
They  are  the  IC,  ICT,  FDTAB,  and  FD  arrays.  The  first  three  array 
names  stored  the  fixed  data  as  discussed  in  the  previous  section. 
These  three  arrays  (IC,  ICT,  and  FDTAB)  are  contained  in  the  labeled 
common  blocks  ARRSET,  TABSET,  and  TABLES,  respectively.  The  informa¬ 
tion  is  available  to  all  programs  in  which  the  common  block  appears. 

The  array  FD  contains  the  input  library  data  used  by  the  program 
to  execute  the  model  for  a  particular  run.  This  array  is  contained 
in  the  common  block  ARRAYS.  The  information  is  passed  through  the 
program  using  the  subroutine  cell  statements.  Integer  values  in  this 
array  are  addressed  by  the  array  ID  which  shares  storage  with  the  FD 
array  through  the  use  of  an  EQUIVALENCE  statement.  In  the  main 
program,  the  data  in  the  FD  array  is  identified  by  pointers  to  the 
section  reserved  for  each  data  type.  These  pointers  are  listed  in 
Table  17. 


TABLE  17.  DATA  IDENTIFICATION  POINTERS 


FORTRAN  SYMBOL 

TYPE  Or  DATA  * 

TYPE  OF  DATA  IDENTIFIED 

no 

ID 

OPERATIONAL  VARIABLES 

IFO 

FD 

OPERATIONAL  VARIABLES 

I  IT 

ZD 

TARGET 

ipt 

FD 

TARGET 

IIB 

ZD 

BACKGROUND 

IFB 

FD 

BACKGROUND 

hot 

ZD 

ENVIRONMENTAL/TERRAIN 

IFEVT 

FD 

ENVIRONMENTAL/TERRAIN 

I  ITS 

ID 

SEARCH 

IFTS 

FD 

SEARCH 

XSETID 

ID 

SENSORS 

ISETFD 

FD 

SENSORS 
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The  data  Is  passed  through  the  call  handler  using  the  following 
procedure: 


CALL  TEST  { . ID(II0),FD(IF0) ,FD(IFT) ,. . . . ) 

and  Is  received  in  the  subroutine  as  follows: 

SUBROUTINE  TEST  ( . IDO,  FDO.FDT, . ) 

DIMENSION  ID0(1),FD0(1),FDT(1) 

Table  18  lists  the  internal  representatives  for  all  the  sections 
of  the  data  used  in  the  subroutines. 


•TABLE  18.  INTERNAL  ARRAY  NAMES 


1  ARRAY  NAME 

TYPE  or  DATA* 

LIBRARY  DATA  ! 

IDO 

IS 

OPERATIONAL  VARIABLES 

ISO 

ID 

OPERATIONAL  VARIABLES 

IDT 

ID 

TARGET 

1ST 

ID 

TARGET 

ZDS 

m 

BACKGROUND 

■SB 

ID 

BACKGROUND 

BEIT 

ID 

ENVIRONMENTAL/TERRAIN 

TDEVT 

ID 

ENVIRONMENT AL/TERKAIN 

ms 

ID 

SEARCH 

IDTS 

ID 

SEARCH 

IDS 

ID 

SENSORS 

IDS 

ID 

SENSORS 

*  ID  -  INTEGER:  FC  -  FLOATING-POINT  NUMBER. 
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IV.  FIXED  DATA  LISTING 


The  following  is  a  listing  of  the  TATAC  fixed  data.  For  an  explana¬ 
tion  of  listing  content,  see  Section  II-D  of  this  Volume. 


IC  ARRAY 


1 

2 

16 

19 

36 

37 

43 

48 

56 

57 

66 

68 

73 

75 

25 

27 

1  1 

t 

113 

136 

138 

161 

163 

1*6 

1  R8 

21  l 

21  1 

21 

1 

21  1 

21  1 

21 1 

211 

211 

21  1 

211 

21 1 

21  l 

21 

1 

21  1 

211 

211 

211 

211 

21 1 

21  1 

21 1 

21  1 

21 

1 

211 

21 1 

211 

XCT  ARRAY 


1 

44 

0 

89 

4 

0 

97 

4 

0 

105 

3 

3  120 

3 

7  151 

2 

4 

145 

4 

9 

214 

4 

9 

263 

4 

7 

302 

4 

7 

341 

6 

25 

522 

6 

2S 

7o: 

4 

0 

71  1 

4 

0  719 

719 

719  719 

719 

719 

719 

719 

FDTAB  ARRAY 


0.00 

.05 

.10 

.15 

.20  .25 

.30  .35 

.40  .45 

,50 

.55 

•  f  0 

.65 

.70  .75 

•80  .85 

.90  .95 

1.00 

1.05 

1.10 

1.15 

1.20  1.25 

1.30  1.35  1 

.40  1.45 

1.50  1-55 

1.60  1.65 

1.70  1.75 

1.80  1.85  1 

.90  1.95 

2.00  2.50 

4.00  5.00 

1.0000000 

.827*345 

.7225450 

.64103*7 

.5742006 

.5177301 

.4691 153 

.4267- 27 

.38936*0 

.3562290 

.3266439 

.3000996 

.2761*39 

.2545597 

.2349471 

.2171109 

.200*517 

.1*599*6 

.  1  *>24041 

.1599404 

.1484955 

.1379713 

.128281 1 

. 11934*1 

.1111041 

.10348*1 

.0964456 

.0899275 

. 3838*99 

.07*2930 

.073100* 

.06*2807 

.063*032 

.0596413 

.0557706 

.0521687 

.0488153 

.0456*15 

. 0427*03 

.0400660 

.0375343 

.0197976 

.00319*2 

.  0009964 

• 

100.  120. 

140.  180. 

.3000  . 

4050  .4675 

.5000 

100.  120. 

150.  180. 

.5130  . 

6330  .6930 

.7000 

5. 

40. 

75. 

1  . 

2.  3. 

-43.  -29.  - 

25.  -39. 

-25.  -20. 

-23.  -17. 

-10. 

.0320000 

.01*0000 

. 07*6000  0 

0000000 

2500000 

.5100000 

2.5400000  6 

.3500000  12 

7000000 

25. 

4000000 

0.0000000 

.000001 4 

.0300042 

0000270 

• 

0000920 

.0002200 

.0005600  0 

.0000000 

00000*6 

0000200 

.0001 400 

.0003500 

.0009000 

0019000 

0. 

0000000 

.0000620 

.0001 

100 

.000*500 

00 1 6000 

0043000 

.0076000 

2. 

>• 

1  • 

2. 

3. 

4. 

1823.  3414. 

4816.  4333. 

366. 

914. 

1829. 

1524. 

1.000 

2.000 

3.000 

4.000 

0.000 

•  t25 

.250  . 

375 

.500 

.625 

.750 

.*75 

1  .000 

0.000 

.250 

.250 

.250  . 

400 

.500 

.500 

.500 

.500 

0.000 

.200 

.200 

.200 

.400  . 

400 

.600 

.600 

.600 

0.000 

.200 

.  .200 

.  400 

.  400 

.600  . 

600 

•  800  1 

•  0O0 

.200 

.200 

.200 

.  400 

.  400 

.600 

. 600*  . 

600 

1.000 

2.000 

3 .000 

4.000 

0.000 

.125 

.250  . 

375 

.500 

.625 

.750 

.*75 

1  .000 

.  1  41 

.061 

.039 

.03* 

02* 

.038 

.053 

.140 

.  462 

.226 

.120 

.057 

.07* 

.050  . 

067 

.091 

.105 

.207 

.154 

.177 

.0*7 

.092 

.071 

.0*5  . 

095 

.  1  49 

.090 

.216 

.113 

.052 

.052 

.046 

.048 

.062  . 

125 

.286 

.2160 

1130 

.0520 

.0520 

.0460  .04*0 

.0620 

* 

«eeeo 


l?SO 
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IC  ARRAY  (Continued) 


1.0000  2.0000  3.0000  4.0000  0.0000  .2500  .5100 
6. 3500  12.7000  25.4000  .<>240  .05*0  .0200  0.0000 
0.0000  0.0000  .9290  .0490  .0210  .0020  0.0000 


0.0000 

.8937 

.0580  .0390  .0080  .0010  .  .0003 

.8290 

.0*40 

.0420  .0050  0.0000  0.0000  O.OOOO 

t. 00000 

2.09000 

3.00000 

4.00000 

0.00000  .25000 

2.54000 

6.35000 

12.70000  25.40000 

.99110  .00790 

.00300 

.00140 

.00030 

0.00000 

.97940  -00790 

.00370 

.00180 

.00040 

o.nonno 

.97710  .00*80 

.00550 

.00320 

.00080 

o.ooono 

.977*2  .00*90 

.00510 

.00290 

.000*8 

0.00000 

1000.00 

1500.00 

5000.00  10000.00 

20000.00  30000.00 

1  .50 

2.00 

3.00 

4.00 

5.00  *.00 

8.00 

9.00 

10.00 

15.00 

20.00  25.00 

35.00 

40.00 

50.00 

*0.00 

70.00  80.00 

100.00 

150.00 

200.00 

.02 

.02  .03 

.08 

.12 

.1* 

.20 

.24  .30 

.73 

1  .35 

2.25 

3.50 

5.40  8.40 

90.00 

90.00 

90.00 

90.00 

90.00  90.00 

0.00 

.02 

.03 

.06 

.09  .13 

.22 

.28 

.34 

.42 

.88  1.65 

5.60 

10.50 

23.00 

80.00 

80.00  80.00 

80.00 

80-00 

80.00 

80. 00 

o.no  o.oo 

0.00 

•  0-00 

.08 

.10 

.13  .16 

.20 

.24 

51  .94 

1.55 

2.50  3.80  5.80 

38. SO  100.00  100. 

00  100.00 

100.00 

100.00  100.00  0.00 

0.00 

0.00  0 

00  '  0.00 

0.00 

0.00  0.00  0.00 

.18 

.31 

49  .72 

1.00 

1.36  2.35  3.85 

10.00 

i6.oo  as 

50  100.00 

ioo.oo 

0.00  0.00  0.00 

0.00 

0.00  0 

00  0.00 

0.00 

0.00  0.00  0.00 

•  20 

.27 

37  .50 

.84 

1.30  1.92  2.75 

5.20 

21.00  100 

00  0.00 

0.00 

o.oo  o.oo  o.no 

0.00 

0-00  0 

00  0.00 

0.00 

0.00  0.00  0.00 

•  27 
26.00 

.34 

56  .90 

1.38 

2.05  2. 88  3.95 

1000.00 

1500-00 

5000.00 

10000.00 

C 

c 

9 

o 

o 

o 

c 

a 

30000. 00 

1.50 

2.00 

3.00 

4.00 

5.00 

6.00 

8-00 

9.00 

10.00 

15.00 

20.00 

25.00 

35-00 

40.00 

50.00 

60.  00 

70.00 

80.00 

100.00 
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100.00 

100.00 

100.00 

100. on 

100.00 

0.00 

0.00 

.19 
12.25 
0.00 
0.00 
7.00 

o.oo 
0.00 
3.90 
0.  00 
0.00 
3.  *5 

95.  120.  150.  1 90.  1.95  1.10  .*5  1.00 
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V.  LIBRARY  DATA  LISTING 

The  following  is  a  listing  of  the  TATAC  library  data.  For  an  explana 
tion  of  listing  content,  see  Section  II-C  of  this  Volume. 
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0.  0* 

10  0  0 

•20S0E+0?  . 1 070E+02 

.8eonE*01  . I600E+00 

»3060E*03 

•  9000E*00  .  moE*o? 

.  JI70E*02  .  U70E+02 

*9000E^O? 

0.  .2500E+02 

0.  o. 

0. 

o.  0. 

0 

0. 8000E-0 1  .3000E+03 

. 9100E*00  0. 

•6000E*01 

0. 

2  2  0  0  0 
• 1 OOOE+OO  . S000E+01 

0.  . 6000E*0l 

•4000E+00 

•  4000£>01  .2500E+00 

0. 

0 

0-  .3600E+05 

. 1 OOOE+03  -2000E+02 

•2550E+03 

•2550E+03  0. 

O 

a 

O 

• 

1  0 

•3600E*05  0. 

0.  .5S00E+C10 

0. 

1  0 

•1333F+01  . 1 860E+01 

.7440i'»0t  -3000E+02 

•3060E-0I 

•I039E-01  0. 

0.  o. 

0. 

1  0 

•lOOOE^OI  . lOOOE+Ot 

. 1 000E*01  .3P0nF*02 

•  4500F  *07 

•2500E+02  .56n0E*01 

. skooe+oi  . ionor*oi 

• 6000E+00 

•4500E-02  .POOOE-P4 

•2000E-OJ  .2000E-08 

•8400E*00 

•  4000E  +  03  • 51 l OE*03 

. S000E+02  .2400E+00 

. 1000E*0I 

0.  0. 

• 

c 

• 

o 

1  0 

• ioooE*oi  . toonr+ot 

.4nnoE*m  . 3000E*02 

•  9250F  *07 

*  1 570E  *02  .S*OOE»Ol 

.  SfiPOE+oi  .!onr»E*r'i 

• 4000E*00 

•4S00E-02  .KOOOE-OS 

.R0n0=:-05  .7000E-08 

.5S0nF*00 

•  aoooF.*03  .snor*o3 
o.  0. 

0.  o. 

o. 

. 1000F40J 

o  o 


FD  ARRAY  (Continued) 
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VI.  SAMPLE  PROBLEMS 


This  section,  containing  sample  problems,  is  provided  to  assist 
the  user  in  setting  up  problems  for  interpreting  the  outputs  of  the 
computer  model.  It  does  not  present  a  comprehensive  list  of  all 
possible  representative  uses  of  the  model.  Each  sample  briefly  des¬ 
cribes  the  problem,  presents  the  execution  list  and  computer  printout, 
and  discusses  the  results. 


A.  VISUAL  OBSERVER  (VI SOB)  -  SAMPLE  PROBLEM  A 

This  sample  problem  exemplifies  the  use  of  the  Visual  Observer 
model.  The  aircraft  search  altitude  is  6,000  feet  with  a  meteorolo¬ 
gical  visibility  of  4  nautical  miles. 

The  specific  variables  which  were  entered  into  the  program  are 
shown  in  the  execution  list  of  Figure  5.  Two  "data  name"  cards  were 
required  with  one  "data  modification"  card  following  each.  All  other 
stored  library  values  remain  unchanged. 

The  results  of  the  computer  run  can  be  seen  from  the  computer 
output  of  Figure  6.  The  cumulative  probability  of  detection  (PAD) 
was  .196  up  to  the  dive  point  and  1.0  at  the  minimum  launch  point. 


S^VS^R  VllVPFp  a  | 

1  =  VISUAL  ODSFFVER 
Pm  F«pHARn-l.«f»*  IMP  INFRARED 
3*  ACT l VF  <  ILi:u*'»MATfP>  TV 
A«  PASS  I VF  tOAYLIGWT)  TV 
5*  FPPVAR0-L'V*K1NG  RAflpp.krri 
6=  FPRPARO-L<lni<I  VP  RADAR*  N3N-MTI 
7»  SYNTHETIC  APERTURE  RADAR 

NAME  I  VAR  VALUE 
PPFR  3  .6000E+04 

ENVI  9  . AOOOE *01 

I  «  MINI PlIM  LAUNCH  POINT 

II  ■  DIVE  BEGINS 

111=  POSITION! 
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VI  ■  CLIPS  TC  ALTITUDE  EEGINS 
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PAD  =  CUMULATIVE  PRORapjlITY  OF  DETECTION 
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Figure  6.  Computer  Output  -  Sample  Problem  A 
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B.  FORWARD-LOOKING  RADAR  (FLR)  -  SAMPLE  PROBLEM  B 


This  sample  problem  exemplifies  the  use  of  the  Forward-Looking 
Radar  model  with  MTI  mode.  The  radar  is  operated  against  a  target 
which  has  a  median  RCS  of  50  M2. 

The  specific  variables  entered  Into  the  program  are  shown  in 
the  execution  list  of  Figure  7.  Only  one  "data  name"  card  and  one 
“data  modification"  card  were  required.  All  other  stored  library 
values  remain  unchanged. 

The  results  of  the  computer  run  can  be  seen  from  the  computer 
output  of  Figure  8.  As  can  be  seen  by  the  flight  profile  printout 
(lines  28-29),  the  aircraft  achieved  search  altitude  at  a  ground 
range  of  22,689  feet  and  the  target  passed  out  of  the  sensor  field 
of  view  at  a  ground  range  of  4768  feet.  The  cumulative  probability 
of  detection  (PAD)  reached  a  maximum  of  .994  which  was  achieved  on 
the  first  "glimpse"  at  22,528  feet.  Closing  to  the  target  did  not 
improve  the  PAD  and  at  a  ground  range  below  8365  feet  the  MTI  could 
not  discriminate  the  target  for  detection.  This  situation  occurred 
due  to  the  following  reasons: 

1.  The  area  of  the  resolution  cell  increases  as  the  range 
decreases,  and 

2.  The  unit  background  RCS  increases  with  increasing  de¬ 
pression  angle  to  target. 
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MAMC  1VAR  VALUE 
TARG  10  .5000P*02 


I  »  MINIMUM  LAUNCH  POINT 

II  ■  01 VF  REPINS 
III*  (DUMMY  ‘'f’SITtONI 

IV  •  TARPFT  PASSES  PUT  0F  FPVfFIXCD  DEPRESSION  AVGLE ) 

V  ■'  SFAFCH  ALTITUDE  ACHIEVED 

VI  ■  CLIP*  TO  ALTITUDE  BEGINS 


X  ■  ALONG *TRACH  GROUND  DISTANCE  TO  TARGET* FT. 
XY  ■  GROUND  RANGE  TP  TARGET  *  FT- 
Tiff-  TIKE  BEFORE  I.AUNCW.SEC. 

PL0S*  PROBABILITY  TARGET  IS  VITWIN  LOS 
PFOV*  PROBABILITY  TARGET  IS  WITHIN  F0V 
P?  ■  SEARCH  TERM  PROBABILITY 
P30  »  DrSCRl«IVAGH.iTY(PETECTI0N) 

P3R  •  DISCPIMINAp  ILITYf  PEOGVITIDNY 

PAD  ■  CUMULATIVE  probability  OF  DFTECTTPn 
PAR  •  CUMULATIVE'  PROBABILITY  OF  RFCOt-N  I T  ION 
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Figure  8.  Computer  Output  -  Sample 
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C.  SYNTHETIC  APERTURE  RADAR  (SAR)  -  SAMPLE  PROBLEM  C 


This  sample  problem  exemplifies  the  use  of  the  Synthetic  Aperture 
Radar  model.  The  search  aircraft  is  flying  at  an  altitude  of  10,000 
feet  and  looking  at  an  offset  of  150,000  feet. 

The  specific  variables  which  were  entered  into  the  program  are 
shown  in  the  execution  list  of  Figure  9.  One  "data  name"  card  and 
two  "data  modification"  cards  were  required.  All  other  library 
values  remain  unchanged. 

The  results  of  the  computer  run  can  be  seen  in  the  computer  out¬ 
put  of  Figure  10.  The  Synthetic  Aperture  Radar  processes  the  doppler 
phase  history  of  the  illuminated  ground  area  over  the  entire  illumina¬ 
tion  period  to  generate  the  display.  All  returns  from  the  target  are 
integrated  into  this  display.  Thus,  only  one  line  of  output  is 
necessary  to  present  the  results.  The  radar  beam  was  directed  26,509 
feet  ahead  of  the  flight  path  and  152,323  feet  from  the  aircraft. 

The  cumulative  probability  of  detection  (PAD)  achieved  was  .178. 
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SENSOR  NUMBER  ■  7 
1*  VISUAL  OPSFRVFP 
?«  FORKApO-LOOHING  INFRARED 
3«  ACTIVE  (ILLUMINATED)  TV 
4*  PASSIVE  (DAYLIGHT)  TV 
5«  FORWARD-LOOKING  RADAR, MTI 
f>  FORWAPO-LOPKING  RADAR. NON-MTI 
7s  SYNTHETIC  APERTURE  RADAR 

NAME  I VAR  VALUE 
(l  PER  S  •  1 500E  +  0A 

0PER  3  . 1000E*05 

X  «  ALONG-TRACK  GROUND  DISTANCE  TO  TARGET.FT. 
XY  ■  GROUND  RANGE  TO  TAPPET »  FT  » 

TIMFs  TIME  BEFORE  LAUNCH, SEC. 

FLOS*  PROBABILITY  TARPFT  IS  WTTWTiN  L«S 
PFOV  =  PROBABILITY  TARPFT  IS  WITHIN  FOV 
P2  s  SEARCH  TERM  PPOPAOILITY 
P3D  «  DISCRIMTNABILITY(DFTFCTIOV) 

P3R  *  DISCS IMINAB ILITYIPFCDPN ITION) 

PAD  «  CUMULATIVE  PROBABILITY  OF  DSTFCTIBN 
PAR  s  CUMULATIVE  PROBABILITY  OF  RFCOGNITION 


X  XY  TIME  PL0S  PFOV  PP  P3D  P3R  PAD 

26509.  J5P323-  .508  1.000  .3SO  1.000  .431  >178 


Figure  10.  Computer  Output  -  Sample  Problem  C 
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VII.  TATAC  PROGRAM  LISTINGS 


This  section  contains  the  listings  of  the  TATAC  computer  programs 
written  in  FORTRAN  IV.  The  following  listing  includes  all  routines 
needed  to  run  the  TATAC  model. 


PROGRAM  TATAC f INPUT. OUTPUT) 

C0MMON/ARRA YS/FD  CP50  > 

C0KKON/ARRSET /  IC  ( SO ) 

C<5»'!KnU/TAPLFS/FDTAP<7?5> 

COPiMCN/T  ASSET /ICTC50) 

cerKAM/RLOCKi /thftah, thftav. phid 
C0KMON/AL0CX3/X,  Y,  X  Y  ,  H,  SR  .  SP.  ANGLO.  AMPLT.  DFLTAT 
C0MMON/PL0CK A/X I . X? . X3 . XA. X5. X6» XVI »XYP»  XY3#  XY4, XY5, XY6 
DIMENSION  IO<PSO> 

EQUIVALENCE  CFO. ID) 

ISNUMSU 
CR=.01 745 
CALL  INPUTI 
Iie  =  ICU 5 
IF0=lCt?) 

I IT= 1C  C3 ) 

1FT=IC<4) 

I IR= IC  < 5 ) 

IFB=IC<6) 

I IEVT *IC  <7) 

IFEVT  =  1C  C  R) 

1 1TS  =  IC  <  9 ) 

IFTSsIC  < 1 0) 

CALL  INPUT2C ISNUKP I 

P0  T0UO.RO. 30,  AO. 50. 60. 705.  ISNUKB 

10  1SFTID*ICCU) 

ISETFOrIC  < 1 2 ) 

ISFNS=t 
ISTYPF  =  | 

P0  TO  000 

80  ISFT 10= IC  < 1 3  > 

ISETFD=ICU4> 

ISF.NS=R 
00  TO  75 

30  ISFTin=lC<15) 

ISF.TFD=1CU6> 

1SF.NS=3 
P0  TP  75 

40  ISFT10= 1C( 17) 

ISFTFD=1CUP) 

JSFNS*4 
C0  TO  75 

SO  1SF.TI0=ICU9) 

ISF.TFD=1C<?0> 

1SFNS  =5 
ISTYPF  *!? 

P0  TO  PO 

60  ISFTin=IC<?l ) 

ISETFD* IC  <PP  > 

ISFNS*6 


VIM 


1STYPF-P 

r-0  TP  *0 
70  JSFTIDMCCPS) 

1SFTFO*  IC  C?a> 

ISFNS*7 

istypf*? 

OP  TP  100 

75  15TVPF»m<!SFTin> 

A5PFCT*FDC  1SFTF0) 

THFTAVeFOf I SFT FO *0 ) 

THFTAH*A5PFCT*Tv'fTaV 
PPTO*FDC  IS^TFOO  5«CP 
OP  TP  POO 

00  THFT*V»FO< ISFTFO* > 0>*C« 

TWFTAP*FOf  1SFTFP*a)*r;p 
PHIO«FD<  1SFTF0*7  >*CP*TKFTAV'5>. 

CALL  FLIOPTC 1SFNS, ISTYPF .FDC  IF0  5 , ID C ISFT in> . FO C 1 SFTFD7 > 
PRINT  POP 
PRINT  POO 
PRINT  POA 
PRINT  PS.XA.XS 
PRINT  PS,XYA.XYS 
R5  FCPKATM  X.3AX.PCF7. 0.1X1) 

PPINT  PO f- 
TlNf*0. 

1PRNTP*  *  1 
I PRV  T  I  *  *  1 
PA  D  *  0  • 

PAR=0. 

90  CALL  PFOPC ISFNS. ISTYPF, FDC 1F0 >» TSTOPi 
1 F  <  I  STOP • FO • I)  f"  T!  POO 
T1PF*TIPF*0FLTAT 
PFPVrJ.O 

call  L"SC ISFVS. ISTYPF, 1  DC  1 1 FVT ) ,FDf IFFVTJ.PLOSi 
PI =PLPS*PFOV 

CALL  SFNFLRC IS* NS, 1r>( I  IT), FDC  TFT)»FDCTFO),FD(  JFFVT), 
t  FDC IFTS I , IOC ISFT ID> ,FDC I SFTFD > , P2 , P3D. P3R 7 
IFCPnn.nT.n, )  r«  to  3*0 
JPENT1 *1PPNT1 ♦! 

!FC< C IPPNTI /? 1+3 l.NF. IPRNT11  00  TO  OO 
PRINT  A?S,X,XY,T INF,PLPS,PF0V,PP,P3D,P3R,PAD,PAR 
00  TO  90 
100  CONTIVUF 
PFOV=t .0 

CALL  SAP c FDC  JFO ) » ID C  1 1 T 5 , FO C I  FT ) , FD C1FP ) , FDC I FF VT > , 

7  ID  C  ISFT  I  D»,  FDC  I  SFT'TI ) .  »’P,  P3D,  P3R  > 

CALL  LOSC  tSFNS. ISTYPF, 1  DC  I  IF VT > , FOC 1 FFVT ) , PLOS > 

PI *PLPS*PFrtV 
PADsPl *PP*P30 
PARePI *P?«P3P 


VI 1-2 


PRINT  203 
PRINT  206 

PRINT  IC5,X,XY«  PL0S  »  PFO V, P2«  P3D,P3R»PAD»PAR 
105  FORMATO  X.2IF7.0,  IX),  6X.7IF5. 3, IX)) 

00  T0  600 

C 

POO  CALL  FLIFHTfISFNS,ISTYPF»FDCIFO)»lDCISFTlD)»FD< ISFTFD) ) 
PPINT  202 

202  FORMAT (|X,/,1X,*I  =  MINIMUM  LAUNCH  POINT*, 

S/,1X,*II  =  OIVF  PFFINS*, 

S/, IX, *111=  (DUMMY  POSITION)*, 

S/,IX,*IV  =  TARFFT  PASSFS  OUT  OF  FOVtFIXFO  DFPRFS&ION*, 

S*  ANC-LF)*, 

S/,1X,*V  =  SFAPCH  ALT ITUDF  ACHIEVED*. 

S/.IX»*Vl  =  CLIME  TO  ALT ITUDF  PFFINS*) 

PRINT  203 
PRINT  204 

PRINT  205«Xl,X2,X3»X/i,X5,X6 
PRINT  205, XY1 , X Y?  »  X Y2 , X Y4, X Y5. X Y6 

203  F0RMAT<1 X,/, 1 X»*X  =  ALONC-TRACK  FP9UND  DISTANCF  TO  *, 

S*TARFET ,FT.+,/,IX,*XY  =  FROUND  RANFF  TO  TARFFT, FT.*, 

S/, 1 X, *T I MF =  TIMF  RFFORF  LAUNCH, SFC.*, 

S/,lX»*PLOS=  FROPAP  ILITY  tarfft  is  within  LOS*, 

S/,1X,*PF0V=  PROPAPILITY  TARFFT  IS  WITHIN  F0V*» 

S/,1X»*P2  =  SFAPCH  TERM  PROPAP ILITY*. 

S/.1X.+P3D  =  DISCRIMI\'APIL1TY<DETFCTI0N)*, 

Z/*  IX,  *P3R  =  OISCPIMINAPIL  1 T Y (RFCOPN 1 T  ION  >*, 

S/,1  X»*PAD  =  CUMULATIVE  PPA»apilITY  OF  OFTFCT I  ON* , 
S/.1X.4PAP  *  CUMULATIVE  PRFBap ILITY  op  RECOGNITION*,/) 

204  FORMAT ( I  X, I 4X,*I*,6X,*I 1*,6X,*I1I*,5X,*IV*, 6X,*V*,7X, *VI * ) 

205  FORMAT  (1X«1CX,6(1X,F7.0)) 

PRINT  206 

206  FORMAT  Cl X»/»4X»*X*»6X»*XY*»5X»  *T IMF*»PX» *PLDS*,2X, *PFOV*, 
S2X.+P2*, 4X»*P3D*,3X,*P3R*,3X,*PA0*,3X,*PAR*> 

IPHASFrO 
T1MF.=0- 
IPRNT 1=1 
X0=O. 

XY0*0. 

TIMEO=0. 

PL©S0=0. 

PF0V0=O. 

ppo=o. 

P3D0=0. 

P3R0=0. 

PADO  =PAO 
PARO=PAR 

CALL  FFOM( ISFNS. ISTYPF »FP(IEO)» I STOP) 

IFC ISTOP.FO. 1 )  FO  TO  600 
TIMF=TIMF,0FLTAT 


VI 1-3 


300  IPHASFelPHASr+J 

|F< I PHASE. FT  •  P  5  G0  T7  600 

IF<  IPHASF.FG.!?  >  PRINT  4!>5, XO » XYO» T tMF0#PL0S0» PF0 VO » 
*  PJ>0#P3O0»P3R0»PADn»PAP0 
IPRNTP*-! 

IFf  ISFNS-FO. 1)  GO  TO  305 

IFOV-IPWAST 

ASPFCT*FD( ISETFD) 

THFTAV*FiKISFTF!>*IFf»V>*CR 
THFTAH=ASPFCT*THFTAV 
ISTYPf=ID< ISETID) 

305  C0NTINUF 

CAUL  SF.ARCHf  ISFMS*  t  STYPE  »  ID(  I  ITI.FDfIFT), 

$  FD<IFFVT)»FDCIFTS)»P?) 

CAUL  F0V(  ISFNS.  ISTYPF,  FDC  IFTS7.PF0V1 

CALL  L0S<  ISEN'S.  ISTYPF,  ID<  I  IEVT5.FDC  IFEVT)>PL0S) 

PI*PL0S*PFAV 

G0  T0 <3 10*3PO>33C'>330*340»340»350>» ISNUPP 
310  CALL  V1S0P  <  FD( 1FT ) #  FfK  1FB) » IPC1 irvT>»FDC I FEVT) » 

.  S  JDf ISFTID>»FD< ISFTFD),P3D>P3R1 
G0  T0  300 

300  CALL  FLTRXFDf I  FT  1  *  FC  t IFP5»FD{ IFEVT7* 

$  IDflSETlDl/FD(ISFTFD>>P3D»P3R> 

G0  T0  300 

330  CALL  TV<  I  SEN'S*  FDC1FT7>FD<IFR)»I0<I  IFVT) »  F0(  1 FEVT  )  * 

*  ID<ISETID)»FP(ISETFD>> I F0  V  #  P3  D  *  P3  R  > 

G0  T0  3P0 

340  CALL  SENFLRC ISFNS,  TOC  I  IT) ,F0< IFT) » FDC IFR 1 . FtK  I FFVT > 
S  FOCIFTS), ID < ISETID:,FDCISETF0>»Pf>»P3D,P3R> 

GO  TO  3P0 
350  CONTINUE 
GO  T0  100 
380  0R  =  P1*PS? 

PAD=P30+CR 

PAR=P3R*0R 

400  !F< ISFNS •  FT*  4 7  F0  TO  410 

1F<<X.LF.X2).AND. C IPRNT! .F0.077  FO  TO  450 
IFKX.C-T.  X?7  •  AND.  CP3D.LT.  .005)  >  FO  TO  600 
GO  TO  400 

410  IF< (PI . OP. PP. OR- P3D > • LT. • 005)  F0  TO  600 
400  IPRNTP= 1 PRNTP ♦ 1 

1F<  <  (IPRNTP/3  >*35  .N'7  .  IPRNT0  1  FO  TO  450 
PRINT  4P5»  X. XY, TINE , PLPS ,  P FOV» PP»  P3P, P3R,  PAD, PAR 
405  F0RPAT < 1 X,0O'7.O» l X) . F5.0»7( 1X»FS.3>> 

IPRNT1 =0 
450  CONTIVUF 
ORN I =0H 
PAON 1 =PAO 
PARN 1 *PAR 

CALL  FFON< ISFNS, lSTYPF.FM IF(1>,  IST0P1 


IFC ISTOP.FO. I I  P0  Tfl  600 
TIMF=TIPF*DFLTAT 
IF  C ISFNS .  PT . 4 1  pa  TO  4*0 

IFCCX-CT.XPI.ANO.  <  IPHASF.FP.  1  ))  PCI  TO  300 
CALL  SFARCWC ISFNS, IS7YPF. ID< I IT>,FDC |FT>, 

S  F0CIFFVT5.F0C IFTSI.P?) 

CALL  F0UC ISFNS, ISTYPE, FPC TFTSI.PFOV) 

460  CALL  LPSC ISFNS, ISTYPF, IDC I IFVTl.FDC IFFVT>»PL0S> 

PI «PL0S*PFOV 

C0  T0C51 0,520, 530, 530, 540.540, SSO), I SNUMB 
510  CALL  VIS3BCF0C IFT), FDC I F8 ) , IDC I IFVT> . FDC IFEVT>. 

S  IDCISFTIOJ.FDC ISETFDl . P3D. P3R> 

G0  T0  SRO 

580  CALL  FL1RCFDC  IFTl.FDC  IFBl.FDC  IFF.VT), 

S  10CISFTin5,FDCISFTF[)),P3F)»P3R> 
t  P0  TO  5R0 

530  CALL  TVC ISFNS.FDC IFT>,FDC IFB), IDC  t IFVT> . FDC  tFF VT >, 

S  IDC ISFTID) , FD< ISFTFDI, I FOV , P3D.P3R > 

C-0  T0  SRO 

540  CALL  SFNFLPCISFNS.IDCIITI.FOC  I  FT  >  ,  FDC  I FB  i  ,  FD  C  I  FF.VT  J  , 
S  FDC IFTS ), IOC ISFTID),FDC I SFTFD > , P2, P3D. P3R> 

C-0  T0  SRO 
550  CONTINUE 

GO  T0  100  •• 

580  PR=PI*P2+C 1 .-PPIkQRNI 
X0=X 
XY0  *XY 
T 1PF0  *T  I  YF 
PL0SO-PLOS 
PF0V0=PF0V 
P?0=P2 
P3D0=P3D 
P3R3=P3R 
DFLTCR=CR-CRN1 
IFCOFLTCR.LF.O. >  PFLT0R=0. 

PAD=PAON I *0FLT0R*P30 
PAR=PARN1 +DFLT0R*P3R 
PAO0=PAD 
PAR0=PAR 
C0  T0  ADO 
600  CONTINUE 
END 

SUBROUTINE  INPUT  1 
C0PNON/ARRAYS/FDCPSO) 

COppon/arrsft/ iccsrn 
COPMON/TAALFS/FOTARC  7?s  ) 

CPt'PON/TARSFT/  ICTC53) 

OIF'FNSION  IDC2507 
FOUIVALFNCF  CFD.ID) 

OATAC ICC  I ), I*1,?4)/I,?,I6, I  9, 36, 37, 43, 4«» 56, 57, 


VII-5 


t  66.6P, 73»75»P5,R7»111»113»136»13R«16|,1 63. 1  R6» I  RP/ 

DATA  (ICO), 1=25, 501/26*21 1/ 

DATA(  ICT( 1 1, J  =  1 » 12)/ 1 . 44, 0»P9,4,0,®7.4,p, 105.3,3/ 

DATAC ICTC I )» 1=1 3.42)/ | 20, 3.7. 151 .2,4, 1 65,4.9, 

S  214’,  4,9,263, 4,7,302,  4, 7, 341 , 6,?5»  522. 6,25. 

S  703.4.0.71 1,4,0/ 

DATAdCTC 15, I =43,50)/? *7  I  9/ 

DATA  10(1), ID(!6), 10(17), 10(1?), !D(36)/0» 10.0.0,0/ 

OAT A  I D ( 43 ) «  10(44).  IDC 45 ), TD(46),  IOC  47 )/?.?. 0,0,0/ 

DATA  IDC56), 10(66), 10(67), 10(73), 10(74) /O, 1,0, 1,0/ 

OATA  ID(R5),I0(R6),ID(111), 10(11?), T0<136)/|,  0,1. 0.1/ 
DATA  10(137), 10(161), I D ( 1 6?), 10(1 R6),ID(lB7)/0» 1,0, 1,0/ 
DATA(FD( I), 1=2, 151/0. ,4000. ,500., 400.  .  400.  ,  400.  , 

$  400. .31 00. . 60000. ,?7., 10., .333,2*0./ 

DATA(FD( I ), 1  =  1 9,35)/?0. 5, 1 0. 7,R. P, .  1 6,306. . .9, 1  I . 7. 

S  1 1 .7,1  1 .7,90. ,0. ,?S. ,5*0./ 

DATA (FD ( I ) , 1=37,4?)/. OP, 300. , . 9 1 , O. , 6. , 0. / 

DATA(FO( I),I=4B,55)/.1,5.,0.,6.».4,4.,.?5»0./ 

DATA(FD( I > ,  1 =5  7 , 65 ) /O . , 36000. , 1 00. , ?0. » P55  •  » 2S5 • , 3*0. / 
DATA(FD( I )» I =6P, 72)/3  6000. , 0. »  0. , . 55,0./ 

DATA(FO( I), 1=75, P4)/ 1 .333, 1 .06. 7. 44,30. , . 0306, 

X  .01039.4*0./ 

DATA(FD( I ) . 1=07,110)/!.. I . , 1 . .30. , . 4597,25. ,5. 6,5.6, 

X  1 . ,. 6, . 45F -P » . 2 F - 4 , . ?E -  1 , • 2E - P , .04,  400.  ,  51  1  .  , 

$  50. ..24,  1  . , 4*0. / 

DATA(FD( I ), 1=1 1 3, 1 OS  1/1 . , 1 . ,4. .30. ,.9?5F7, 1 5-7,5. 6, 5. 6, 

S  I . ,.6». 4SF-2, . OE-F, . OF -5, - 7F-8, . 55. 600.  ,  511  .  ,0.  , 

X  0.,t.,3*0./ 

DATA  (  Ft)  (  I  ),  1  =  130,1  601/2  4.  ,  1000.  ,  .  ?E- 1 2,  3  600.  ,  60000.  ,  .03, 
S  90 • , 1 0 « «?.5,9P.,3P.,.2F-5».5«.P,0.»0.,0.,0.,.9S, 

S  0..15..0. ,0./ 

DATA(FD( 1 ), 1=1 63, 1051/24. , 1000. , . 2E- 12, 3 600. , 60000. ,. 03, 
S  90 .»10.,2.5,90.,30.,.2F-5,.5,.0,0.,0.,0.,0.,.95, 

X  0*,15.,0.,0./ 

DATA(FD(  1  ),  1  =  1  00,2  01/6.  .  1  000.  ,  .  ?E-  12, 3  600.  ,  60000.  »  50.  , 

5  15.,0.,.03,0.,0.,30.,.1--5,.5».0,0.,1.»10.,.95, 

$  36. ,36. ,36 . , 1 . / 

OATA ( FO ( 11,1=211 ,250)/ 40*0. / 

DATA(FOTAR ( I ), 1=1 .44)/ 

*  0.00.0. 05,0. I O, 0. ■ 5,0.20. 0.25,0. 30.0.35. 0. 40,0. 45, 

X  0.50. 0.55,0. 60, 0.65, 0.70,0. 75,0.R0,0.R5,0.9O,O.°5, 

X  1.00,1.05,1.10.1.15.1  20,1.25.1.30,1.35,1.40,1.45, 

X  I • 50. 1 • 55  » I • 60, 1.65,1. 70,  1.75,  1  . RO, 1 . RS.  1 . 90,  1 . 95. 

X  2.00,2.50.4.00,5.00/ 

DATA(FDTAP< I), 1=45. PR)/ 

X  1 .0000000, 0.5279345, O. 7225450, 0.641 03R7.0. 5742006, 

S  0.5177301,0.4601 153, 0.4? 67127, 0.3903690, 0.3 5 62290, 

X  0.3266439,0.3000996,0.2761 930, 0.°5 46597, 0.2 3 40471, 

$  0. 2171109. 0.200 9617,0. 1P599R6, 0.1 7? 4041, 0.1 6 90404, 

X  0.I4R4 95 5,0. 1379713,0.1 ?9PR II. 0.11 03 4P1, 0.1111041, 


VII-6 


t  0*  I034RRI  ,0.09644S5,0.0399P7S,0.0938R99,0.0782930, 

$  0. 073 1  008. 0.0632807 ,0.0639032, 0.059641 3, O.OS577 06, 

S  0.05?  1  687,0. 04  9. 91S3»0.<V  56915,0.  0427803, 0.  0400660, 

S  0.037  53/13.0. 01  97976.0.0C3  1  9RP.  0.  000996  4/ 

DATAFFOTAP  (  I ).  I  =  P9, 104  >/ I  00. , 1  00.  .  I  AO.  .  |  ftp.  ..3.  .405. 

S  .4675. .5. 1 00.. 120., 150., 130., . 5 1 3 , . 633 . . 693 » • 7/ 
DATA<FDTARC  I  >.1  =  105.11  9>/S. .  40.. -75.  >  I  .  ,2.  .3.  . -43.  .  -29.  , 

S  -25  ..,-3  9.  ,-25..  -20-.-23.  .-1  7..-10./ 

DATAFFDTAPC  I >» 1=120, 150>/.032. .01 R, .0086,0.  ,  .25. .51 . 

S  2-54,6.35,12.70,25.40,0., 1 . 4F- 6, 4. 2F-6, 2. 7F-5, 9 . 2F-5, 

S  2.9F-4,5.6F-4,n. ,P. 6F-6.P.0F-5. 1 .4F-4.3.5F-4.9.0F-4.  • 

S  1.9F-3.0..6.PF-5. 1 . 1F-4.R.5F-4, 1 . 6E-3, 4. 3F-3 . 7 . 6F-3/ 
OATACFOTAPF I>» 1=151, 1 64>/2. , I . 3. .4. . 1823. .341 4.. 

S  481 6. ,4333. ,366. ,914., 18?9. , 15?4./ 

OATA<  FDTAF!  ( I >, 1  =  165, PI 3>/l .3..4..0.. . 1 ?5» .25, . 375, . 5. 
$  .625, .75, . R7S  » 1 . , O. , 3* . 25 , .4, 4* . 5 , 0  • , 3*. P, 24. 4. 34. 6, 

S  0.»2*.2,2*.4»2*.6,.R,!.  ,  0.  »3*.P,2*.4,3*.  6/ 

DATA  <FI>T AR  c  I  > ,  1=21 4,P62>/| . , P .  , 3 .  , 4.  , 0 .  » .125,  .25,  .375,  .5, 
S  . 62 5 , .75, •  875  » 1 ., . 1 41  . . 061 ,.039, • 038, .0?P, «03R, .053, 

$  .140, .462, .226. . 1 PO, . 057, . 07R, . OSO, . 067, .091 , .105, .207, 

S  .  1  54,  •  1  77,  .097,  .  092,  .071,  .035,  .095,.  1  49,  .  09,  . ?.  1  6,  .113, 

S  .052,  -05P,  .04  6,  .  04R,  •  06P,  .125,  .21* 6/ 

DATA<FDTAP< I >, 1=263,301 >/1 . ,2. ,3. ,4.,0. . .25,-51, 

S  2.54, 6 • 3S  » 12. 70,25. 40, . 924, .056,-02,4^0.  ,  .929, .04P, 

S  . 02 1.- 002. 3*0.,. R937,. OSR,. 039, .OOS.. 001,. 0003,0.,. 889, 

$  .064, .042, .005, 3*0./ 

DATA  <  FDTAF!  (I),  1=302, 340>/l. , 2. , 3.. 4. ,0.,. 25, .51, 

S  2. 54, 6. 35. 12. 70, 25. 4,. 9R11 , . 007R, . 0064, . 0030, . 001 4, 

S  .0003,0.  , .97  94, .0079, .0063, . 0037, . 001 P, . 0004, 0. ,. 977 1 , 

S  . 0063, . 0066, . 0055 , . 003? , . 0003,0. , - 97762, . 0D69, . 0068. 

S  .0051, .0029, .00063,0./ 

DATAFFDTARC I >, 1=341 , 430 >/ 1 000 . , 1 500. , 5000. ,10000., 

S  20000 .,30000.  , l ., 1 *5,?. ,3. ,4. ,5* ,6., 7., 9. ,9., 10*, 

$  I  5. ,20* ,25. , 30. , 35. , 40. , 50. , 60. , 70. ,90. ,90. , 1 00. , 

S  150., 200., .02, .02, *03,. 06,. OR, *12, .16, .20, .24,. 30, 

S  .35, .73, 1-35. 2. 25, 3. 5, 5. 4, 9. 4, 21. 5, 90. , 6*90. , 
t  0* ,.02, .03, .06,. 09, . 13,. 17,.P?,.?R,.34,.^J,,pp, 1,65, 

S  3., 5.6, 1 0.5,23. ,90. ,7*«0. ,5*0. , .08,. 1  ,.  13,.  1  6/ 

DATA' FDT  A3 < 1 >,  1  =  431,521  > / . 20, . P 4 , . 5 1 , . 94, 1.55,2.5, 

S  3.8, 5.3, 13*, 33.5,  100., 5 *100.,  10*0., .09, .19, .31, .49, 

*  .72, I. ,1*36, 2. 35, 3. 35, 6. 3, 10. ,16. .25. 5, 1 OO. , | 00. , 

S  12*0. ,.14, .20, .27, -37, .50, .34, 1.3,1.92,2.75,3.75, 

S  5. 2, 21. ,100., 14*0. ,.22, .27, .34, .56, .90, 1*38, 2.05, 

$  2.88,3. 95, 12. ,26./ 

DATA CF0TA9 f I ) ,  1  =  522,61 1  >/1000. , 1500^,5000.  ,10000. , 

S  20000., 30000., 1 ., 1*5, 2., 3., 4., 5. ,6., 7., 3., 9., 10., 
f  15. ,20. ,25 . ,30. ,35. , 40. , 50. , 60* ,70. ,30. ,90. , 1 00. , 

J  150. ,200. ,.01 , .03,. 04,. 09,-1 6,. 24,. 33*.  43,.  53,.  65, 

$  .78,1 .7,3* I , 5 • 2, 9. P, IP. 25, I  7. 75,38. ,34.  , 1 00. ,5*100., 

S  0. , . 02 , .04, .03, . 1 5, .22, .32, . 44, .56, .7, .38,2,07, 4. 1 , 
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f  7.4, 12.5,20. 5.32. »  72.  ,  1  00.  »  6*  1  Of).  ,  5*0.  ,  .  07,  .  11  *  .  1  5# 

*  .19/ 

DATACFOTAP( I >. I  =  M  2. 702  5/.?$, .32..R, 1 . A,?. 25. 4. 2, 

S  7«6,t?«25»30.»f6.,1 00. » 5+100. *  1 0+0. ,.1,.?1,.4,<67, 

5  1.05. 1.55. 2. 23. 4. 1,7. .11. 2. 1R. 2. 31. 5.  AO. ,120. .120.. 

S  12*0.  ..  14,  .22.  .34.  .  4ft  .  .  ^4.  I  •  1  .  |  ■  75.  2  •  A  ,  3  •  *7.5. 4. 7. 6. 

S  34. . I  AO* . I  4*0. . .2?* •  ?  R, .3ft. . 5B. • P9. 1 .3. 1 . 9.2. 7.3. 65. 

$  15. .47./ 

DATA(FDTAP(  I  ).  1=703 .710/95..  120.  » I5O..1R0.  » 

S  1 .95,1 . 1 . . 65. 1  •/ 

DATA(FDTAR(  15.1=7  11 , 7 1 R5/9S. , 120. , 1  SO.  , 1  AO.  . 

S  2. 7 A, 1 . 72. 1.1.!./ 

DATAfFDTAP  <  15.  1=71  9. 7255/7*0./ 

RETURN 

END 

SUBROUTINE  INPUT? < I SNUMP 5 
C0MW5N/ARRAYS/FDC2SO5 
C0MMON/ARRSET/ !C(5P5 

01MFNSION  10(2505. XTEMP( 105, 1 TEMPI (10), TTFMP2 (10,105, 

■ #TEMP3< 10,105 

ECUIVALENCE  (FD.IP5 

DATA  XEXFC , XFVDR, XPPER. XTARG/4HFXEC. *HFNDR»  4H0PFP. 4HTARC/ 
DATA  XBACK.XFNVI,XSFAR,XSENS/4PPACK,4HFNVI,4HSEAR,.4HSENS/ 
DATA  XPLANK, XTERR/ 4H  , 4HTFRR/ 

ISNUMR=1 
NTEMP 1 =0 
DO  I  Js 1,10 
1  XTEMPC.J5=XRLANK 
READ  5.AAA.KSNUWP 
5  F0RKAT(A4, AX, 1 1 5 

IF (AAA • EO • XEXFC  5  00  TO  15 
PRINT  10 

10  FORMAT < *  EXECUTE  CARD  MISSING  OR  OUT  OF  9RDFR* 5 
00  TO  30 

15  1 F  < (KSNUMP . GT • 0  5  «  AND . ( KSNUMO . LT • B  >  5  00  T0  25 
PRINT  20.KSNUMB 

20  FORMAT ( *  INVALID  SENSOR  NUMBER--*, 125 
00  TO  30 

25  ISNUMT?  =  KSNDMR 
30  1SFT  =  5+ I SNUMP 
35  READ  40.XNAMF, IN 
40  FPRMAT(A4,2AX. 125 
11  =  0 

IF (XNAME .F0.XFNDR5  00  TO  200 
IF (XNAME • EC • XOPER  5  11  =  1 
IF(XNAMF.EP.XTARG5  11=2 
IF(XNAMF.FO. XPACK5  11=3 
IF (XNAMF . FO • XFNV I  5  11=4 
1F(XNAmF.F0.X~ERR5  I  I  =  4 
1F(XNAME.EC.X1'.FAR5  I  1=5 
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IFtXNAMF.FO. XSFNS)  II=ISFT 
IF<  1 1  • NF • 0 )  CP  Tf I  SO 
PRINT  45.XNAMF 

45  FORMAT (*  INVALID  DATA  NAME--*»A4> 

C0  TP  35 

50  ILPOPsO 

NTFMP 1 =NTF  MP 1 ♦ I 
XTFMPCNTFMPl  ) s XNAMF 
ITFMPI (NTFMP1 )  =  IN 
NTFMP? =0 
11*2*11-1 
T2=2* 1 1 
13=?* 1 1  + 1 
III=IC( II )-l 
I  I2  =  IC  < I?)-l 
II3=IC<13) 

60  IL0OP  =  ILO3P+ 1 

IF< ILOPP.PT •  IN  }  G0  T0  35 
RFAD  70. 1VAR, VALUE 

70  F0RMAT(RX,  IP.F10.4) 

IF  < I VAR . GT • 0 )  C0  TO  80 
PRINT  7S. I VAR 

75  FORMAT (*  VARIABLE  NUMBER  0UTSIDE  RANEE  0F  DATA  NAME* 
$*  ARRAY--*. 12) 

E0  T0  60 

BO  IF  ( I VAR . LE #(113-11. ?))  G0  T0  100 
PRINT  90. I VAR 

90  FORMAT (♦  VARIABLE  NUMBER  0 UTS  IDF  RANEE  0F  DATA  NAME* 
S*  ARRAY--*. 12) 

E0  T0  60 

100  J  J=  I 1 1  +  I VAR 

NTEMP? -NTFMP? ♦ 1 
ITEMP2(NTFMP|  ,NTFM3?5  =  JVAR 
TFMP3  ( NTFMP I »  NTF MR? ) = VALUE 
IFUVAR.CT.  (  I I2-II 1  ))  GO  T0  110 
ID<JJ)  =  IFI X( VALUE ) 

GO  TO  60 

110  FD(JJ)=VALUE 
GO  TO  60 

200  PRINT  205. ISNUMR 

205  FORMAT( I  X./, IOX,*SFNS0R  NUMBER  *  *»I1, 

S/.IX,*  1=  VISUAL  ORSFRVFR*. 

S/.IX,*  2=  FORWARD-lOOKING  INFPARFD*. 

S/.IX.*  3=  ACT  I VF  (  ILI-UMINATFD)  TV*. 

S/.IX,*  A-  PASS  I  VF  (DAYLIGHT)  TV*. 

S/.IX,*  5=  FORWARD-LOOKING  RADAR, MTI*. 

S/.IX,*  6=  FORWARD- LOOK  INC  RADAR ,N«N-MT I *, 

S/.IX,*  7=  SYNTHFTIC  APFRTURF.  RADAR*) 

L00P*0 

210  L00P=LOOP*  1 


IF  (LOOP,  GT. 1 0 )  CO  TO  240. 

IFC  XTFMP  CI.',I*P  )  .  \  F  «  XBLANK  >  CO  Tfl  215 
RETURN 

215  IFCLPOP-FO. I >  PPINT  220 

220  FORMAT  f  1  X*/,1X,«NAMF  IVAR  VALUE*) 

NN  =  ITFMP  l.CLO^P  ) 

00  230  J=1.NN 

PRINT  225* XTFMPf LOOP)* 1TEMP2CL00P* J)»TFMP3  CLO0P* J> 

825  FORMATCI X»44,3X, I2*3X*E10.4) 

230  CONTINUE 
P0  TO  210 
240  PRINT  245 

245  FORMAT  <  *  MORE  MODIFICATIONS  THAN  LISTFD*  > 

RETURN 

END 

SUBROUTINE  FLIGHT! ISFNS, ISTYPE*FDO» IDS.FDS) 
COMM0N/PLPCK1 /THFTAH,THFTAV*PM1D 
C0MMON/PLOC><2/'UL*HP*SPD»SPL*SPC»SPP 
C0MM-0N/BLOCK3/X,  Y» XY. H, SR* SP» ANGLO » AVCLT* OELTAT 
C0MM.ON/PLOCK4/X1  » X2 , X3 » X4, XS* X6* XY I . XY2* XY3 * XY4,  XY5*XY6 
DIMENSION  FD0 (20 ) *  I DS C 20 > * FPS C20 ) 

C8«. 01745 
Y«=ED0C1) 

ML«FD<K2) 

SPL=FD?(S>*1  .  AP7778 
!F<ISFNS.PT. 4)  GO  TO  500 
MC-LD=FOV(iO>*Cf* 

SPD<=FD3<4)*1 .6P7778*CPSCAVCLD> 

SRI =  FD0  <  R  > 

OELTAT  =  Fnr  f l p  > 

XY2«HL/TA\ (ANGLO) 

X2«S0RTC  XY2*XY?-Y*Y) 

TEMPt*ATANCY/X2) 

XY!=SRt*CCS<AVGLD) 

X1=XY1*OOS(TFMP1 ) 

IFCISFNS. FO • 1 )  TO  TO  300 
I F  <  I  ST  YPr  •  EO  »  ?  )  C-0  TO  200 
100  PH1D=AN0U> 

HP=FD0C3 ) 

5PP=FD0 { 7  >  *  1 . AR777R 

ANGLC  =  F00f  1 1 )*CP 

SPC  =  FD0(6)*1 .AF777R*C0S(  ANGLC ) 

SR2  =  FPP(<?) 

XY3«XY2 
X3  =  X2 
XY4*XY2 
X4*X? 

XY5=SPF!TCSRP*SP2-WL*HI.) 

X5  =  SCRT<  X  YSiXYS- Y* Y) 

X6»XS*CHL-HP>/TAN(ANCLC) 
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XY6«S0RT!X6*X6fY*Y> 

GO  T«  400 

200  CALL  COVFR!  1SFNS,  IDS#FDS.  Y»HL»XMAX»X*«IN> 
1F!XKAX.PE.X? >  P0  T">  210 
PH 10= ! AT AN ! HL/XP )  *TWFTAV/2« >/CR 
PRINT  20S.PHID 

205  FORMAT ! *  DEPRESSION  ANCLF  TOO  STFFP.  PHIO  RESET' 
S  F5.2,*  DFPREFS  «  > 

PHID=PHID*CR 

CALL  COVER! ISFNS* IDS, FDS,  Y.  WL.XMAX#  XHI*t> 

210  X5=XMAX 

XY5=S0RT! X5«X5+Y*Y> 

SR2=SCRT!HL*HL*Y*Y+XY5*XY5> 

X6  =  X5 
XY6=XYS 
X3  =  X? 

XY3=XY2 

'  X4=AMAX1 !X2»XMIN> 

XY4=S0RT!X4*X4*Y+Y> 

HP=HL 

5PP=SPL 

SPC=SPL 

ANGLC=0. 

GO  TO  400 

300  CALL  COVER! ISENS, IDS.FDS. Y,HL#XPAX,XMIN> 

GO  TO  210 
400  X=X1 
XY*XY1 
SR=SR1 

ANGLO *AC0S!X2/XY2> 

ANGLT=ANCLO 

H=XY*TAN!ANOLT) 

SP=SPD 

RETURN 

500  IF! ISFNS.E0.7)  FO  TO  600 

CALL  COVER! ISFNS# I DS»  FDS#  Y#HL»  XMAX»XWIN> 

X4=XMIN 

XY4=S0RT!X4*X4+Y*Y> 

X5=XMAX 

XY5=S0RT!X5*X5+Y*Y> 

GO  TO  700 
600  CONTINUE 
700PX=X4 

XY=XY4 

H=Hi. 

ANGLO =ATAN!Y/X4> 

ANGLT=ATAN!H/XY> 

SR =5 CRT !XY*XY*H*H> 

SP«SPL 

DELTAT«F0O ! 1 P ) 
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RrTUPN 
t  NO 

SUBROUTINE  COVFRC ISFNS* !DS#FDS#Y»H.XMAX#XMIN> 

FXTFRNAL  FUN1 

C0MN0N/PL0CK1  /TVFTAH.  THFTAV*PHID 
C0MM0N/BLOCK5/X2YAX,XIMIN,VISMAX 
DIMENSION  1DSC20>.FDSC20> 

ce**ot74S 

IF<ISFNS.FO.  1  >  CO  TP  ?00 
XMAX*H/TAM  CPHI0-THFTAV/2. 1 
XMIN*H/TAN  CPH10*THFTAV/2. ) 

TEMPI  «TAN  C THE T AH/2 •  > 

YMIN»XMIN*TFMP  1 

IFCY.PT. YMIV)  PO  TP  t?0 
RFTURN 

ISO  YMAX«XMAX*TFMP1 

IFCY.PT. YMAX1  PO  TP  150 
XMIN*XMAX“ C YMA  X- Y  J/TFMP 1 
RETURN 
150  Y*0. 

PRINT  155.Y 

155  F0RMATC*  0FFSFTCY)  CREATFR  THAN  Y*AX.  SET  Y=  *»F10.4> 

RFTURN 

800  VISMAX«FDS C I  1 
RMAX=VISMAX 
ICMASK«IDSC1  1 
IF<YI  210.220.230 
810  YrAFSCYl 
PRINT  2 l 5 

815  FORMAT C *  ONLY  POSITIVE  OFFSETS  C0NS1DFRFD.  USED  ABS.  VALUF* ) 
GO  T0  230 
820  XMAX*RMAX 

TEMPI *FUN 1  CO. 5*CP 
xmin=h/tanctf*pi > 

IFC  ICMAStC.  EG  •  0. 1  XMIN=0. 

GO  TO  300 

830  IFCY-LT.RMAX5  PO  T0  250 
XKAX=0- 
XMIN*0. 

GO  TO  300 

850  XMAX=SCRTCRMAX*RMAX-Y*Y> 

XMIN=0. 

IFC1CMASK.FG.0)  GO  T0  300 
GUESS =AS  IN  C  Y/RMAX  )/F.  P 
TEMPI =  H/TAN  C  FUN  I <PUFSS)*C8) 

TEKP2*Y/S INC  PUFS5*CB ) 

IFCTFMP1 .|_E.TFM^2>  GO  T0  260 
XMAX*0. 

XMIN«0. 

GO  TO  300 


260  GUFSS*rifFSS*5. 

IFCfUFSS.GT. oo.  )  P!IFSS«oo. 

TFPP1  *W/TOM(  FI  INI  <r.HFSS»*C«> 

TFKPP»Y/siN<ri)Fss*r«> 

IFCTEPP1 .GT.TFPP? >  GO  T0  P70 
IF< GOFSS.LT. 90*  I  GO  TO  260 
XP1N=0.  ‘ 

GO  TO  300 

870  CALL  APPROXCfUESS.TWETAD, IFAIL) 

‘  TFMPI*FUN1  (THFTADXCB 
TFPP1 =H/TAN<TFPPI  ) 

XP1N=S0RT<TFpP1*TFPPI -Y*Y> 

300  X?PAX=XPAX 
XI  P IN*  XM IN 
RETUPN 
END 

5UFR0UT INF  APPROX (  GUFSS, ANS  »  I  FAIL! 

FXTFRNAL  FUMF,FUMDF 

1FAIL*0 

ANS*GUFSS 

00  100  1*1,50 

TFPP1 *ANS 

AMS  *  AM  S  -  (  FUN'F  (  AN.S  I  /  FlINOF  C  ANS  >  > 

TEPP?=A° SCANS -TF-V!P1 > 

IF CTFPP2.LE. I . E- A I  RETURN 
100  CONTINUE 
IFA IL* I 
RETURN 
END 

FUNCTION  FUNF(THFTAD) 

EXTFRNAL  FUN  I 

C0PP-0.V/PLeCKP<'HL,HP,SPD,SPL,SPC,SPP 

C0PP.0N/PLOCK3/X,  Y,XY,H,  SR,  SP»  ANGLO.  ANGLT.DFLTAT 

CR= • 01 745 

THFTAR=THFTAD*CB 

TEPP l =FUN 1 <  THFTAOI *C8 

FUNF=HL/TANCTFPP I >  -  Y/S INCTHFTARJ 

END 

FUNCTION  FUNOFCTHFTADI 
EXTERNAL  FUN  I , FLNP 

CBPPON/RLOCKP/HL , HP. SFO, SPL.SPC.SPP 
COPPON/PLOCK3/X, Y» XY,H, SR, SP, ANGLO, ANGLT.OELTAT 
CP* . 01 7  AS 
THFTAR=THFTAD*CP 
TFMP1 *FUNI <THFTAn)*CP 

FUNf)F  =  -HL*FUNJ>(THFTAD>*CR/'SINCTEPP!  >**2 
S  ♦HL*COSCTHFTAR)/S IN  CTHFTAR  >**{» 

FNO 

FUNCTION  FUN  I <  DECREE > 

FUN  I « - . 003 1  A? *DF GREF *DEGREE*. 4?309*DEGREE* 1 2. 9176 


FND 

FUNCTION  FUNP(f)F.PRFF  > 

FTNP«-.OOAPR4*OEGPFP+.  4P309 

END 

SUBROUTINE  GEOMf ISFNS.  ISTYPF.FD'''.  ISTflP) 
C0MM0N/PL0OXP/WL.  HP.  SPO>  SPL»  SPC  ,  SPP 
C0MM0N/PL0Oi<3/X»  Y»XY.H»SRi  SP» ANGLO*  ANGLT.  DFLTAT 
C0MM0N/BL3CK4/X1  *X?*X3*X4,X5#X6*XYI.*XY2*XY3»XY4*XY5>XY6 
DIMENSION  FDO  < 20 ) 

IST0P*0 
CB». 01745 

IF < ISFNS.PT . 4 )  PO  TO  POO 
ANGLD=FDO<  1  0  5*CR 
IFtXY.PT.XYP5  PO  TO  130 
XY«XY»SPO+OFLTAT 
TFMP1*ATAN< Y/XP5 

X*XY*COS  <TFvp ] 5 

H?XY*TAN (ANFLT  5 

IFt XY*LE • XYP  5  P0  TO  IPO 

TC'0  =  tXY-XYP5/SPD 

X*X?*SPL*TG0 

SP*SPL 

H=HL 

110  XY*S0RT tX*X+Y*Y5 
ANGL0«ATAN< Y/X5 
ANC-LT=ATAN<H/XY5 
IPO  SR*S0RT<XY*XY*H*H5 
RETURN 

130  IFttlSFVS.FO.  1 >. OR.  GISTYPE. EO-255  GO  T0  POO 
ANGLC«FD3  < 6  5*CR 
IFtX.C-T.XS)  GO  T0  140 
X«X*SPL*DELTAT 
IF<X.LF.X55  00  TO  110 
TP0x t X -X5  5  /  SPL 
X*X5+SPC  +  TC-0 
SPrSPC 

HxHL-SPC  *T GO *T4.N  (  AVPLC  5 
lFtX.LF.X£5  00  TO  110 
TGO* t X-XA  5 /SPC 
X«X6*SPP*TP0 
SP=SPP 
H»HP 

P<3  TO  110 

140  1FCX.CT.XM  GO  TO  150 
X«X+SPC*DFLTAT 
SP=SPC 

H»H-SPC*OFLTAT*TAN ( ANGLO  5 
IFtX.LF. XA5  GO  TO  I  I  0 
TG0«tX-XA5/SFO 
X«X4*SPP*TG0 
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SP=SPP 

H«HP 

00  TO  MO 
150  XsX*SPP*DFLTAT 
G0  TO  1 1 0 

200  X®X*SPL*DFLTAT 

IF(X.LF.XS)  GO  TO  t  1 0 

IST0P=1 

GO  TO  110 

FND 

SUPROUTINF  SFARCHG ISFNS. ISTYPF, IDT.FDT# 

SFOFVT.FOTS.PP> 

c«) KMC M/ FLOCK  1  /THFTAH.THFTAV,  PHID 
COMMON/PLOCK3/X,  Y.XY.H.SR.SP.  ANGLO » ANGLT, DFLTAT 
C0MM0N/BLOCK4/X! » X?, X3, XA, X5. X A, XY1 , XY?» XY3, XYA.XYS. XYG 
C0MM0N/PL9CKS/XPMAX.X1  MI\>,  VMS* AX 
COMK0N/PLOCK6/ALPHAH, ALPHAL, ALPHAV 
DIMENSION  IDT(20),  FDT  C  PO  )  .  FDEVT  (20)  .FDTSC20) 

F(D1»D2)=P.*ATAN((-01-1. +  SORT  ( (D1 *1 • >**2*Dl*D2*D2>  >/(01 *02 ) > 

CPI=3.M1S9 

DTX  =  FOT  < 1 > 

DTY=FOT  (?) 

DTZ=FDT(3) 

NT  =  IDT  ( 1 > 

TSPA.CFsFDT  (10) 

DSX=FDTS(2> 

I»SY=FDTS  f  3  ) 

ANPTC=CPI/2. -ANGLT 
TANGTC  =  TAN(ANC-TC) 

IFCAVCLO.LF.  (ATA,VCDTY/DTX>>>  G0  T0  10 
TEMP 1=S IN (ANGLO) 

OTXXrDTY/TEPPl 
C-0  TO  20 

10  TEMPI =C?S(AN0L9) 

0TXX=0TX/TFMP1 
20  TEMPI =TA"PTC*TANGTC 

TEMP?=(DTXX+r)TZ/TAN(ANGLT>  )XH 
ALPHA V=F (TEMPI »TFMP2) 

ALPHA  H=2  .*ATAN((DTX*SIN(ANGL0  ) ♦DTY*CSS (ANGLO >  >/(?•* SR >) 
ALPHAL  =  ?.*ATAV((DTX+COS(A.VCLP  >*SIN(ANCLT>* 

$  OT  Y*S  IN  ( ANGLO  )«S  IN  (ANGLT  )*r)TZ*COS (ANGLT  >  >/(2 . *SR>  > 

If(X.LT.XA)  GO  TO  P50 
AT*DTX*DTY*SIN(ANGLT>  *■ 

S  (OTX+DTZ*SIN(A\iOL0)  +  DTYM)TZ*C0S(ANCLO)>*C0S(ANGLT> 
AT=AT-»(NT*(NT- I  >*TSPACF/0TX  > 

IF( ISFNS.FC. 1 )  GO  TO  100 
IF< ISTYPF . FO . 1  )  PH  I D  =  ANGLT 
PHIC=CPI/P. -PHID 
TPH1C  *TAN ( PP 10  > 

IFOISX.LF.O. >  00  TO  30 
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tempi «TPwrr*rpMrc 

TEMPPeDSX/H 
8ETAV»E<TEMP  1 » TFMP? ) 

IFtBETAV.LT.THETAV)  G0  TO  40 
30  BETAV=THETAV  , 

40  IFtDSY.LE.P. >  GO  TPI  SO 

BETAH*P.*ATAV<OSY/(P.*SR»i 
IFfBETAH.LT.THETAH)  GO  TO  60 
50  PETAH=THFTAH 

60  0MEGA=PETAH*<S!NCPHID*RFTAV/8» >-S INtPHID-BETAV/? .  >  ) 
SRPH10  =  H/SI'V<PHID) 

AS*0PEr-A*SPPHlD*SFPHID 

00  TO  POO 
too  P'LOC^FPTS  f  4) 

TEMP  I  eXPNAX  +  SOPT  (  X?MAX*X?MAX+  Y*  Y*H*H> 

TEKP2»X1  MIN+SCRT  f  XI  M  IN*Xi  M I N* Y* Y*H*H1 
AS*H*»’L3C*AL0f  (TFMPI/TEMPP) 

800  ARATIO=AT/AS 
G*F0EVT<6J 

•  TEMPl*t700./G)+ARATI3*DELTAT 
IFtTEMPl .LT. IEO. >  OC  TO  270 
850  PP=1.0 
RETURN 

870  P2*l . -EXPt-TEMPI ) 

RETURN 

END 

SUBROUTINE  LOS ( t SEN* , I  ST YPF  » IDEVT ,  FDEVT.  PLOS  > 
C0MM9N/TAPLFS/ EOTAF C  725 ) 

C0MM0N/TAPSET/ ICT  C  50 ) 

C0MM3WPLOCX3/X,  Y,XY»H»SP,SP.AMGLO.ANGLT.DEUTAT 

C0MM0N/RL0CK4/X1  ,  XP,  X3.X4,x5>X6,XY1  . XY? » X Y3 , XY4 , XYS, X Y6 

DIMENSION  I DFVT  C  ?0  > .»  FDF  VT  C  PO  7 

I E  <  X  •  LT •  X 4  )  GO  TO  35 

C8=.0t  745 

ICFL0S=I0EVTC3> 

ANGLM=EDEVTC2 )*CA 
RAT10M=FDEVTC:n 
PNMA$K=I .0 

IF(AS1N<H/SR7.GT.ATAN(RATI0M7>  GO  TO  10 
PLOS=0. 

RETURN 

10  1FCANGLM.LF.0. 7  GO  TO  ?0 

PNMASK=PNMASX*C I . -FXPt - AN0LT/ANCLM7 ) 

80  PCELOS=  I  .0 

IFt 1SENS . OT «  4 )  GO  TO  30 
IFtTCELOS.EO.rn  GO  TO  30 
ZENITH=A.NGLT/C:P  +  90. 

JJ= ICT  <4  7 
NX= ICT <5  7 
IX«JJ 
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IY*IX*NX 

IF<H.LF.5000.  )  F9I  TO  25 
JJ= ICT <7  ) 

NX=ICT<BI 
IX*  JJ 
IY*IX*NX 

25  CALL  INTFRI  CFOTABC  IX>..FDTAPCIY>»NX»ZFNITH»PCFL0S> 

30  PLOS=PNMASK*pr.FL0S 
RETURN 
35  PL0S= 1 • 0 
RETURN 
END 

SUBROUTINE  INTFR1  (XX.  YY»NN,X,Y> 

DIMENSION  YYM  ).XX(|  > 

Y=yr<i ) 

IFCNN.FO* 1 >  RETURN 
DO  I  K=1»NN 

IF6X.LE.XXOO)  GO  TO  2 
]  CONTINUE 
Y*yy<NNi 
RETURN 

2  IFCK.FO.n  RETURN 

Y= < <  X-XX { K- 1 nX<XXOO-XXCK-l > n*CYYCK>-YYCK-I ))*YYCK-1) 

RETURN 

END 

SUBROUTINE  FOVC ISFNS, I STYPF, FDTS, PFOVI 
C0MM0N/BLOCK1/THETAH,THFTAV,PWID 
C0MM0N/RLOCK3/X, Y, XY.W,  SR , S P» ANGLO , ANPLT* DFLTAT 
COMMON/RLPCKA/XI  , X? , X3.X4. X5»  X6* XYI »XY?»XY3»  XY4»  XYS»  XY6 
DIMENSION  FDTS  <20) 

CPI *3- 1 41 59 
S I GX=FOTS (51 
S  I G Y=FDTS  <  A 1 

IFdSFNS.FO.  1  )  GO  TO  100 
IF< ISTYPE.E0.2)  GO  TO  100 
IE<X«LT • X4 )  GO  TO  100 
ANGTC=CPI/2. -AMGLT 
IF<SIGX.GT.O. )  GO  TO  40 
PFOVXrl.O 
GO  TO  50 

40  XV  =2 .  +SR*TA\' <THFTAH/?»  >/ < 2 . B2B42B*S I GX> 

PFOVXsFRRFUNCXV) 

50  IFXSIGY.GT-O. )  GO  TO  60 
PFPVY= 1 .0 
GO  TO  90 

60  I FC <ANGTC*THFTAV/2 • > • LT • (CPI/2. )>  GO  TO  70 
PY1 *.5 
GO  TO  FO 

70  XU  =M*<TAN<ANGTC  +  THFTAV/?. :-TAN<ANGTC»»/(2.B2R4?P*SIGY> 
PYI *• 5+ERRFUN ( XL  1 ) 
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*0  XLP»R*f TAN  <ANf;T0'l-TAN<  ANRTC-THFTAV/P.  >>/<?.  P2«a2R>*S  ICY) 
PYP».5*FPPFUN<XLP> 

PF0VY«PY1  *PY2 
90  PF0V»PFOVX*PF<?VY 
RETURN 
100  PFAV*|.0 
RETURN 
END 

FUNCTION  FRRPUN < X > 

C*l. 1PR379I 6709551 
ERRFU-N*- t • P50 
IFCX)  6.1,1 
1  X?*X*X 

1FCXP.LT. 1  PO- >  C0  TO  10 
DENS«0. 

C-0  T0  1  1 

10  DENS=C*FXP< -X?) 

11  IFCX-2.25 )  A, 4,2 

C  C0PPUTE  ERRFUN  US  IMP  A  CONTINUED  FRACTION  EXPANSION 

?  N=76./X-4. 

A I *XP*  C  XP*4. 5 ) ♦?. 

AP*X*CXP*CX2*7. >+R.?5> 

Bl*X*CX2*CX2+5. >*3.75> 

B?=X?A  C X?  +  C  X?  +  7 .57+11 . P5>+1 ♦ P7S 
T*3. 

□0  3  1*3, N 
A3*X*A?+T*A1 
A1  =A2 
AP=A3 

B3*X*B?+T+Bl 
B1*B2 
B2=B3 
3  T=T  +  •  S 

CERRF*DENS*A3/<P.*P3) 

ERRFUNrl . -CPRRF 
RETURN 

C  C0PPUTE  ERRFUN  PROP  THF.  MXLAUR 1 N  POUFR  SERIES 
A  N= 1 1 • 5+X+5 . 

T*5. 

ERRFUN  =  X«  < I . -X2*C3. 33333 333333333F- t -XP* 

*  C. 1 -XPACP. 3POo5?3P095?38E-2-X2*A. 62962962962963F-3  )  >  7  > 

Y*X*X2*X2*X2*X°/P 1 6. 

D0  S  1=?»N 

Y»-XP*CP. *T -  1 . >*Y/<T*<?.*T*1 . >7 
ERRFUN=FRRFUN+Y 

5  T*T+  1 

FRRFlINrC+FRRFllN 
CFRRE  = 1 • -ERRFUN 

6  RFTURN 
FND 
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SUBROUTINE  ATMOS <FOT,Fr>n,  JDFVTjFDFVT,  I!>S»FDS,  ISFNS* 
STAT  #CA  »CM) 

C0MM0N/TARLFS/FDTARC7P5  >. 

common/tarset/ictcso) 

C0MMON/RLOCKI /THFTAH.THFTAV. PHID 
C0MM0N/nLOCK3/X>  Y»XY>H» SR. SP. ANGLO. AA'PLT* PFLTAT 
D 1 MEN S  ION  FDT  <  PO  > , FDR  f  PO  > , I DF  VT <PO >  »  FDFVT  ( ?0 >  , 
SIDSfPO>,FDSC?0> 

VG=FDFVTt4) 

IFtlSE.VS.F 0.?>  FO  TO  POO 
RT=FDTt4> 

RB  =  FOB (  1  5 

CI=ABStRT-RR)/AMAX! (R7.RB) 

TEMP  1 = • 00 1  *W 

TEMPI =.44/ (TEMPI  *1 .?>* . 56*EXPC-.08*TEMP1  >  +  .073 
VS=VF/TEMPI 

IFC1SENS.F0.1  >  XLAMDA=FDSC4> 

IFt ISFNS.F0.3)  XLAMDA=FDSC I5> 

IFCISF\'S.F0.4)  XLAMD4--FDSC15) 

S I  CMA=  t  3.9t?/tVS*A07F.  )  >*  CXLAM.DA/ .  55  > 
if<isfns:fo.  t  )  GO  TO  5 
GAMMAT  =  FDS  <  9  ) 

IFtISFNS.FO.O)  GO  TO  100 
5  TAT=EXP<-SIGMA*SR> 

I TAT  =  1DFVT  <4> 

IFC 1TAT.F0. 1 >  GO  TO  300 
10  ICLOUDrlDFVTtt ) 

RTE=FDFVT(! ) 

TEMPI = . ? 

IFtICLOUD.FO.?)  TFMPlr.6 
IFt 1 CLOUD. FO.  3  >  TEMPI *1. 

ZK=TEMP 1 /RTE 

TEMPI  =  1 . *ZK*<1  .  -TAT) /TAT 
50  CA=C I / TEMP  1 
55  IF (  ISFNS . FO. 1 )  RETURN 
60  CM=CA/<?. -Cl > 

CM=1 . -{ l . -CM)**C AMMAT 
RETURN 

JOO  GA IN  =  FDS  C 1 9 ) 

1ST  YPF  =  1  OS  1  1  ) 

1F< ISTYPF . FO. 1  )  PH  I D= ANGLT 
TEMP  1  =  PM I D -THFTAV/P. 

TFMPP-PH  >*TPFTAV/5>. 

SRMlNrH/SINf-TF^Pl  ) 

SRMAX-H/S  iN'fTFMPP) 

ZMIN=?.*S!  GMA  *S!’M  1 N 
ZMAX  =? . *  S I G MA*SRMAX 
TAT=FXP( -ZMAX) 

JJ  = ICT  < 1  ) 

NX* ICT  <? ) 
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X  Y®  I X+NX 

CALL  INTFR  1  < FDTAP (  I  X  7  *  FOTAO  <  I YJ  »  MX.  ZMIN  *  FOT^-IN) 

CALL  INTFP1  <FDTAQC  IX>.FDTAPCIY)»NX.ZMAX.F?Zi*AX> 

XNA  =  (  CA  J\*S  ICMA*SPMAX+.?5  )* 

S  <SRMAX/SRMIN*F?ZMIN-F?ZMAX) 

XNT=XNA*RT*TAT 
XNB®XNA  +  RF-*TAT 
•CA=ABSCXNT/XNP-|  .  1 
pei  Te  to 
POO  RH=FDFVT (  5  ) 

DELTB=FDR<2  ) 

TEMPI®. 00013*H 

TEMPI =< I  ./TFMPI  5*<  1  . -EXPC -TEMPI  5) 

TEMP2=.0002R-*H 

TEMP2=< 1 • /TFMP2 ) *  C I .-FXPt-TFMPP) ) 

XV0°I  •  A3?*RP'»'‘XPt  .  0652* C  OFLTP-273.  7  > 

TEMP2=<  .  27P*TFVP27/  (  VC*6.  07*  >  +  .01  7*XV!B*TFMP1 

TAT*FXP( -TFMP2*  SP  +  .001  7 

RETURN 

300  I AZIM* IDFVT  <57 
HH=H 

SR1000=.001»SR 
8MAX®AMA  XI (  F T » PB  7 
ZENITH=ANOLT/.Ot  7/5+90. 

IE<IAZIM.F0-P7  GO  T0  350 
JJ® ICT <317 
•  NX® ICT  <323 
NY® ICT  <33  7 
IX=JJ 
IY=1X+NX 
I Z= I Y+N Y 

CALL  INTFR2  < HH /SRI POO. FDTAP 1 1 X ) . FDTAB ( I Y 7 , FDTAP (  I  Z  ) 
S  NX.NY.DPP. XXXX.NX7 
JJ=  ICT (37) 

NX® ICT (38) 

IX=JJ 

IY=IX+NX 

CALL  INTER  I < FDTAP < I X 7 , FDTAP ( I Y7 , NX. ZFN I TH, RF7 
00  TO  360 
350  JJ= ICT  <3  A  7 
NX  =  ICT (35  > 

IX=JJ 

IY=IX*NX 

?Z=IY+NY 

CALL  INTFRP :hu,SR | 000. FDTAP ( I  X 7, FDTAP C I Y 7. FDTAP C IZ) 
$  NX.NY.OPR.XXXX.NO 
JJ=ICTCAp7 
NX® ICT  <  A 1  7 
IX=JJ 
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IY*IX*NX 

CALL  INTFR1 (FOTAPC 1 X 3 . FntAP C I Y3 »NX. ZFN1TH,RF3 
360  CA=CI/C1  .♦nPR/<RMA»RF3  3 
00  T0  55 
END 

SUBROUTINE  V ISOR ( FDT*  FDE, T  DFVT#  FDFVT* 

SIDS»FPS.P3f>>P3R3 

C0MM0N/BLOCK3/X,  Y»>;Y»H,SR,SP»ANOL0»ANPLT»DFLTAT 
C0MM0N/nLOCK4/X!  ,X8,X3»X.4.  XS,  X6» XY1  .  XYS>» XY3  »  XY4,  XY5»  XY6 
COMM0N/BL0CK6/ALPHAH»ALPIJAL.  ALPHAV 
DIMENSION  FDT(203>FDBC?03» I DFVT  (P03»  FDFVTCP03* 
SIDSCP03.FDSCP03 
ISFNS  =  1 
THFTA=.R 
CR  =  «  01 7  45 
XK1 *.0265 
XK2=.P4 
XK3=.44 
XK4=l-6 

CT*XKI *THFTA**XK?+  t  XK3*TPFTA*+XK43/ 

S  (  (AMIN  1  (ALPHAH.  ALPHA V 3+60. /CP 3**P 3 

call  ATM0S(FDT,FDB, IDFVT,FDFVT,XXX,FDS*ISENS, 

SXXX.CA#  XXX 3 
TEMPI =CA/CT 

1 F(TEMP1 . LT • • 05  3  00  TO  90 
I F (TEMPI . LT> I «  3  03  TO  40 
AKr.43 
00  T0  50 
40  AK=- .57 

50  TEMPI *4. 2+ ( (TEMPI  - I • 3**P 3 

P3D= .57+AK+SGRT < l . -EXPC-TFMP1 3  3 
A=3  . 

B=2.9F-6 

1 F( X . OF.  X?  3  00  T0  55 
VALPHArO. 

00  TO  60 

55  VALPHA=(SP*SORT( H+H+ Y*Y3/(SR+SR3 3/CR 
60  TVA* ( A+B* ( VALPHA  +  +  3  3  3/P • 

TFMP1 = ( (  AMIN  1 (ALPHAH. A LPHAV3+60. /CR  3/TVA  3**2 
IF(TFMP1 .LF.3.P3  0-0  TO  100 
TFMPl  =  { (TEMPI -3.P  3**P  3/ 1  l . 

IF(TFMPI .LT.  I  RO.  3  00  TO  70 

P3R=1 .0 

RETURN 

70  P3R*1  . -F.XP(-TFMPI  3 
RETURN 
90  P3D=0. 

100  P3R=0. 

RETURN 

END 
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SUPPmjTINF  FL!R<FOT<FDn.FnFVT< I OS  *  F0S<P3D<  P3R  J 
COMMON/Pl.OCKt /TPFTAW,  TPFTAV.  PH  ID 
C0MMON/RLOCX3/X.  Y<XY<H<SP< SP<ANCLO<ANGLT<DFLTAT 
C0MM0N/RL0CK6/ALPHAH, ALPHAL< ALPHAV 
DIMENSION  F0T(?0)»r DP<20J<FDFVTC20J< 
SIDS<FOJ<FDS<FOJ 
ISENS=2 
DFLTB=FDBC2J 
EMM8=FDR<3 J 
DELTT=FPT (5  J 
EMMT  =FOT  <  6  > 

G=FDEVT  <  A  J 
ZK1=FDS<SJ 
ZK2=FDS<6 J 

CAUL  ATMOS  <  XXN<FDB<XXX<FDFVT<XXX<  XXX< ISENS< 
STAT<XXX.XXXJ 

XNV=AMIN1  <ALPHaH,ALPKALJ/AMAX!  <ALPHAHtAUPPAUJ 
IX*1 

too  TEMP ! =0. 

XN*2.*THETAV/AMINI  calphah,alphavj 
IFCXN.C-T.SOO.  >  G0  T0  320 
IF<IX<NF<2J  02  TO  HO 
XN= A. *XN 

lFCXN.OT.SOO. J  GO  TO  320 
GO  TO  130 

HO  IF<XN<GT<300<J  GO  TO  120 
SNRT=6< 1R*FXP<-. 00252* XNJ 
GO  TO  ISO 

120  S'NRT=?<  901  76-.  OOP* < XN-300.  > 

GO  TO  ISO 

130  IF<XN<GT<300<J  GO  TO  1  AO 
SNRTcR  <?*F.  XPC*.002  AP*XN  J 
GO  T0  150 

1  AO  SNRT=3<R967l - . 002+ < XN -300. J 
150  SNRT=SNRT* < 1 < ♦ . OOP*0*G > 

IF<XN.LT< 100. J  Go  TO  160 
SNMRT  =3. - <  002*XN 
GO  TO  170 

1 60  SNMRT =5. -<o??*xn 
170  XMiRT  =  ZX1  +FXPC  ZK2*XN  J 

TDIFF=APS  <DPLTT'*FPMT  -DFLTU+FMMP  J 

SNRD=<SNMRT*TDIFF'  TAT/XMRTJ*t  CDFLTP*FM.MB/300<  1**3  J* 
$  SDRTCl </<7.*XNVJ) 

SNR=SNRD-SNRT 
IF<SNR<LT<S<J  GO  TO  310 
TEMPI =1 .0 
GO  TO  320 

310  TFMP 1 = 1 < -FPRFNN<A0S<SNR>/1. A1 API AJ 
TFM.P1  s.5*TFMf>| 

IF<SNR<GT<0<J  TEMPI *1. -TEMPI 
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320  IFUX.FP.23  GO  TO  390 
330  P3D=TFMP I 

IFITEMPl  .LF.O.  >  09  T0  390 
IX  =2 

GO  TO  tOO 
390  P3R=TFMP1 
RETURN 
FND 

SUBROUTINE  TV< TSFNS#  FDT»  FOP# IDFVT#FDEVT, 
SIDS^FOS, IF0V,P3Di P3R3 
COMMON/PLOCKl  /TWFTAH,  THFTA V*  PB  JO 
C0MMAN/ALPCK3/X. Y,XY»U»SP»SP.ANGL9»AN0t_T»DFLTAT 
C0MM0N/PL0CK6/ALPHAH, ALPHAL. ALPHAV 
DIMENSION  FDTC203.FDBC20), I DE VTC203 #FDFVT C20 3 # 
SIDS<?0t»FDS<20>  ' 

ASPECT  *FDS  < 1  3 
BAND=FDS  C  5 1 
DIAO=FDSf 65 
PNUM=FOS (6+1  FA V  5 
GAMMAT *FDS  C  9  3 
XFC=FDSC1 1 5 
XI=FDSC1P3 
xiMAX=rns<i35 
XIP=FDSCIA) 

Xl_AMDA  =  FDS ( 1 S  3  ‘ 

XNL=FDS (165 
XNR=FDS  <17  5 
6T=FDS(20) 

RTE=FDEVT  < 1 3 
G=FDFVT<63 

CALL  ATMOS  <  FOT»FOO  » IDFVT*  FDFVT » IDS»  FDS» I SENS* 
STAT»XXX,CM3 

EARFA.  =  AS  PEG.  T*D  I  AP*D  J  AG/  (  1 . +ASPECT*ASPFCT 3 
RESP0N  =  1  0.-k*(CAL0(-10(XI  :  - ALO C- 1  0 ( FARFA 3  3 /GAMMAT 
S  -AL9G10<.67A*XFC3 3 

IFLAG=0 

IF<ISFNS.F0.33  GO  T0  10 
ICLOUDr i  r VT  C  t 3 
I SUN AN  = 1 1  -  VT  <  2  3 
TRANSM.=  EDS<  103 
CFAC= I . 

1F< ICL0UD.E0.23  C-AC=3.\635 
IF(ICL0UD.FO.3>  C*AC=tO. 

HS=595 • 

1F<  1SU.NAN.F0.23  HS=465. 

IFOSUNAN.FO.33  H  9=PAS. 

XICPN=RFSP0N*HS+RTF*TRANSM*TRANSM/(4.»CFAC3 
GO  TO  20 
10  POUT  =  FDS  < 1 P  3 

ASTR*THFTAV*THFTAH 
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20  IF <  I5FNS .  fp. 3  >  r-n  to  30 

XIAVF=FAREA*  (  X  TOON/  f  FMIJK*FNUM>  >*PAMMAT 
PCI  TO  AO 

30  XICONP  =  RFSPPN  +  PTF*P0llT*TAT*TRANSM/C4.i<ASTR*SR«'SR> 
XIAVF=EAI?FA*  f  X1C0N?/  CFNIJMtFNUM)  >**GAMMAT 
40  IFCX1AVF.LF. C.«S*X IMAX>1  GO  TO  50 
FNUK=l.At4P1  4*FMIJM 
IFLAG=1 
GO  TO  00 

50  IFCIFLAP.FO.O)  GO  JO  A 0  . 

PRINT  55  » F\'tJM 

55  FORMAT  <  *  TUPF  S  ATt  RAT  ION  *  NFV  FNIIM*  *,F5.?> 

40  XVO ISF=XIP+XIP/CP. *B AND  > 

XF=SORTCXVL*XNR/  C  1  .  A  t  4?  1  A* ASPECT )  > 

XN0E=1 .?4*XF 
XNFT  *«3 I *XF 

XNEL»S4S  .  *D  I  AG/  C  XLAMDA*FNUM*SORTC  1 .*ASPFCT* ASPECT  >  ) 
XN0L=3.6A*XNFL 

XNV*AMIN  1  <A|_PHAW>  ALPHAL  J/AMAX1  CALPHAH, ALPHAL > 

IX*1 

100  TEMPI =0. 

XN*?.+THFTAV/AM!NI  C ALPHAH. ALPHAV) 

IFCXN.GT.700. >  GO  TO  3P0 
IFCIX.NE.P)  GO  TO  110 
XN=4.*XN 

IFCXN. GT • 700* )  GO  TO  320 
GO  TO  130 

110  IFCXN. GT. 300. 7  GO  TO  IPO 
SNRT=A. 1 P*FXPC-.00P5?+XN> 

GO  TO  150 

120  SNRT  =2*90176-.  OOP*  C  XN-300. 7 
GO  TO  150 

130  IFC  XN. GT .300. 7  GO  TO  1  AO 
SNRT=R.2*FXP<-.00248*XN> 

GO  TO  150 

140  SNRT=3.P9f.71-.r>0?*CXN-300.  7 
150  SNRT  =  SNRT*  < I .+.003*G*G7 
ZNFT  =  C XN/ < XV V  +  XNFT  7  )**P 
ZNEL=CXN/CXNV+XNFL7  7**2 
PS  I Y=SORT < 1 . ♦Z\FT«ZNFL> 

GAMMAY=PSI V/SPRTC 1  .  *7.N FL*P.  *ZNET7 
IFCXN. GT. CXNOE/3.  7  7  GO  TO  P10 
BE TAT  =  C . 3*XF /XN 7  *ERRFUN  C  3.  *XN/XE 7 
00  TO  PPO 

210  BETAT=CXF/XN7*C .0RG*.?33*ERRFUNC2. 45*XN/XF77 
220  IFCXN.LF • tO. 7  GO  TO  PAD 
1FCXN.0T.XN0F7  GO  TO  P50 

'  le-1 
RSOF=0. 

2P5  I»I*P 
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IFfl.GT.SO)  G 0  TO  ?AO 
XIN*I*XN 

l  F<  X  IN*  GT •  XN51F  >  GO  TO  PAO 
PHIX  =  XIN/XN!1L 

RLN=. A37*(AC0S<PHIX>-PHIX*SCPTC1 *-PKIX*PHIX>> 
IF<XIN.GT* <XN0E/3. 5 >  GO  TO  POO 

RSOF=RSOF*f.R!  I  *RLN*£XP< -4. 4*XIN*XIN/<XF*XF  >  )  >/t  I*  I  ) 
GO  TO  RP5 

830  RSGF=RS0F*<.637*RLH*F.XP<-3-*XlN+XlN/(XE*XF.)))S(l*l> 

GO  TO  P35 
840  RSOF-1 .0 
GO  TO  PAO 
850  RSOF-O. 

GO  TO  320 

960  SNRD=(?.*RSOF*CP*GT*XIAVE/XN)*SORT<. l*XNV/< ASPECT* 
s  PS  I  Y*<OT*CT*BFTAT*GAMMAY*XIAVF*|.  AF.-19*XN01SE>  >> 

SNR=SNRD-SNRT 
IF<SNR*LT*5. 1  GO  TO  310 
TEMPI *1-0 
GO  TO  320 

310  TFMP1  =  1 • -FRPFUNt ABS  (SNR>/1*414P|4> 

TEMP1=.S*TEMP1 

IFfSNR.GT.O. )  TEMPI =t. -TEMPI 
3P0  I F (  I X . F.O •  ?  1  GO  TO  350 

IF<RS0F.GT*0.  >  GO  TO  330 
P3D-0. 

P3R=0. 

RETURN 

330  P3D=TEMP1 
1X=2 

GO  TO  100 

350  JFCRSOF.GT.O. >  GO  TO  390 
P3R=0. 

RETURN 

390  P3R=TEMP1 
RETURN 
ENO 

SUBROUTINE  SENFLRC  TSENS* IDT*  FDT * FDB*FDEVT» 

S  FDTS » 1DS»FDS*P?*P3P*P3R> 

COMMON/TAPLFS/FDTABf 7P5) 

C0MM0N/TABSFT/1CTC5O) 

CPMMON/PLOCK 1 /THFT  AH*  THFTAV*  PHID 
COMMON/PLOCKO/X*  Y,  XY » H » SR * SP* ANGLO ,ANGLT» DFLTAT 
DIMENSION  IDTC201*  FDTCPO) , FDP 120> * FDEVTC20) » 

S  FDTS<20)*1DSC20!*FDS(25) 


CPI =3* 141 59PA5 
C2«.304R 
C3=l • 6R9 


r>  n 


CB».0I745 

SMORbO. 

VC«0. 

NT* IDT (1  ) 

PN»D. 

•  VRTbO. 

•  VB»0. 

c 

G*FDFVTC6> 

VTT*FDT  1 127 
C«  3.F  8 
H«H*C2 

VG=SP/t . 687778 
G0*FDSC2> 

PRF«FDS<4> 

PX*FDS (5  > 

FLAKDA  *FDS ( 6 ) 

WS'FDS ( 7 ) 

PHIM=PHID-TBFTAV/2. 

WbFOFVT  C  7  > 

TAUPV.'=FDS<I2> 

TAOPl  =F[)E  (13) 

'TAUP2=FDSf I  4) 

VF*FDS  CIS) 

VPRFsFOS (16) 

CAS«FDS<I2> 

ISEASbIDSCI > 

DTXbFCTCI >*C2 
DTY*FPTC2>+C2 
PS JbFDT ( I  I >*C8 
XXXXbO. 

VP*FDECS> 

GG  *  0. 

PNR  =  FDS (3  > 

C  SELECT  PA\D 
1PAND* 1 

IF<FLA VOA . LT • C  .0?5>> IPAVD=2 
IF<FIAI«'QA.LT.(.0I25>>  IBANDbD 
BAND«  IPA.VD 
C  CGVFRAOF 

TbT*C2 

C 

P3DtO. 

x«x*c? 


GFO-TTPY 

SR*SP»C? 

PHIbA'ITLT 
ocy=sr*tmftam 
DGXbC*TAUF'V’/{2.*C0S(PHJ  )  ) 
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TSPACF=FOTCtO>*C? 

RNJC-NT 

opki*j*amimi  cdgy,ooy> 

RNJC  l*00f',IN/{NT*C  TSPACF»DTX>  > 
ratsipbamivm  cnN.ic.nNJCi ) 
SIGFACbAKAXI  (  1  .  .RATSIG) 

XTErAFINI  (CTX.DGX) 

YTE*APIN1  <CTY.DGY> 

AC—DCXTDCY 
ATE'XTF*  YTF 

ABE*AKAX1 {  <A0-ATF),0.  > 
TEPPl=C',lS<PHn**? 

TEK.P  1 =SORT  <  SR*SR*TFI^P1  -Y*Y> 
ALPHA=ATANP< Y.TFPP1 ) 
VRT*VTT*COSCPHn*COS<PSI-ALPHA> 
TFFP1 =S INCPHI)**? 

TEMPP  =  SINCPHlM>*-*5> 
TEFP3*SCRT<C3S<PNI  >> 
TFKP4=SORTCCOSCPHIfO  ) 
PAF=r.0*TF^PP*TEPP3/(TFMPI*TEKP4> 
C  RESOLUTION  GEOMETRY 
AG1=AG 
ABE 1*0. 


IF-<ABE.EO.O.  >CQ  T0  POO 
100  0X=AMAX1  (DGX.DC-Y5 
DN=APIN 1 <DGX»DGY> 

1F<0X.LE.I'P)  GO  TO  POO 
IFCDiV.LF.V'P)  GO  T0  110 
AP«V*P*(DX*DM)/B. 

G0  T0  IPO 

110  AP*DN*  <OX*t‘'P  5  /P  • 

IPO  AL*AFAX1 C (AG-AP>,0. ) 

FI *1 • 0 
TEFPI =Ft 
F1*F1/TEPP1 
A1 *AP*Ft  +AL 

1F< ATF . GT . AP>  GO  TO  130 
ABE1 *Al -ATF 
GO  TO  140 
130  ATFL=ATF-AP 

ABF1 =A1 -FI  * ATEL-AP 
1  AO  AG1*F|*AC- 
C  RADAR  CROSS  SECTION 
POO  TFFP1 sSINCPHI  > 

ATrDTX*DTY 
l  =  A+IAANf) 

SIGTFrFOTf I ) *f . 3/FLAMDA  >**. 3 
SIGTF»1O.*AL0G1O(SICTF)*PHI/(CB*1S.» 
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TFHPS'PSI 

ir<PSl. PT..7R5i»)  TFPPS'CPl/P.  -PS  I 
S K-TE'S  I CTF -S/4S.  *PS  I/CP 
SIGTF'10. *+<STPTF/l 0. ) 

SIPTF'SIGFAf>SlC-TF 

eio  tepp?=ag/at 

IF <TFPP2  •  CF  •  l  •  >  pn  T!1  ??0 
S I  PTE'S  I GTE *TE  PP2 

P20  PHI'PHI/CS 
jj'iCTao) 

NX'ICTO  1  > 

N  Y'  ICT  <  1 2  > 

IX'JJ 

IY*IX+NX 

IZ'IY+NY 

CALL  INTER? ( PHI , PANO, FDTAP C IX) , FDTAPt I Y  V* FDTAP ( I Z > , 
S  NX»NY»SIC-A|  »  XXXXi  VX 1 

siPni*io.**fsips! /in. > 

882  PHI 'PH I *CS 

SICBF'TEPPl  * ( APF 1  *SIFRt  } 

S  IGP'TF.KP  1*(AG1*F.  1PB1  > 

CALL  ATTEN ( IBAND> CAS>  H* ISEAS. PHI , TAUAR > 

TAUAR* 1 . / ( 1 0. ** (TAUAR/ 1 0. > ) 

310  PNB'O. 

RECEIVER  POVfR 

TEMP3'4*CPI*CPI*>SR»SR 

PDfl'TAUPl *PX*FAF*TAUAR+CPI/TEKP3 

PRTT'S I PTF  *PDO 

PRBT'S ICPF+PDn 

PRRB'SIPP*Pnn 

TEMP 1  =TAUP?  *PAF *TAUARhFLANDA**P/  C  A.  *TFHP3  ) 

PATT'TEHPI *PRTT 

PART =TFYPl*PRPT 

PARB'TFMPI *PRPn 

TFKP1 =THFTAH*prf/(WS*CR> 

TFPP1 'SORT (TFNPl  ) 

PTT'TFfPI *FATT 
PTB'TFtfPl  *FAPT 
PRB'TEMP I *F  ARB 
SNRAT  =(PTT*PTR>/ CPRB+PNR) 

1F< ISENS •  FC • 6 }  GO  T0  361 

CLUTTER  FREQUENCY 

VR'O. 

330  TEPP I 'S INI  PH  I > 
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TFMP?=cns<PHi) 

TFKP3*SIN(ALPHA) 

TFMPAsO^St ALPHA) 

VVY=C9*03  *V0*TAUPl-'*TFMP4+TFMPl**3*C/(TFMP?*FLAMr>A*H) 

VVX=C?*C3*VG+THFTAH*TFMP3/FLAMDC 

Vk‘=2*  *Vk'/FLA.vl>4 

VS=FOSf 7)*CF/(S0PT(2.*CP1 )*THETAH) 

VC*SOPT  C  VVY+ *?*VVX**?*VV  **2*VR**?*VS**2) 

CLUTTFR  ATTFNIJATION 

VT*SRRT(Vr.**?*VF**P  +  VPPF**P) 

VTH=F0S<P3) 

VB=FLAMOA*PBF/(P. *03*0?) 

IF  C (2 .  *VTH  )  •  LT •  VR  )  f-01  T0  333 
PRINT  33P.X.Y 

332  FOPMATf  I  X.35HTHRFSH3LD  SPFFD  IS  GREATER  THAN  1/2. 

S12H  RLINO  SPEFD»5X,2HX=E12.5,?X»2HY=F12.5> 

P3D=0. 

G0  T0  370 

333  FN1=(VRT-VTH)/V2 
N1*FN1*. 00001 

IF(FN 1 . LT •  0*  )  N 1 =N 1  - 1 

N 1 =N 1*1 

FN1=N1 

TFMP1=(VRT*VTH)/VR 

IF(FN1 .GF. TEMPI )  C-0  TO  334 

P30=0. 

GO  TO  370 

334  .  FTH=2. *C2*C3*VTH/FLAMDA 
•  FMF  =  FDSC1  ) 

IF  (FMF)  9,7.9 
7  FMF  =12. 

PRINT  F 

R  FORMAT  <5X.  42HAN  MF  VALl'F  WAS  NOT  INPUT  SO  MF  IS  SFT«12.) 
9  F 1 1  =  FMF/6. 

FKK  =  |. /Cl. 414  *FTH**FII) 

FRT  =  P.+C2  *C3  *VF  T  /FLAMDA 
FBK  =  FTH  *1.414  **< 1 • /FI  I ) 

GG  =  FKK  **?  *FRT  **C?.*FII) 

IF  <FRT  .GT.  FBK)  GG  =  I. 

PTT1  =  PTT  *GG 
TEMPI =0« 

CALL  GAIN (TEMPI , FBK, TEMP?, VT.O.XXXXJ 
TEMP?=.5-(TFMP?/(VT*S0RTCP.*CPI )>> 

IF  (TFMP2  .LT.  0.)  TEMP?  =  0. 

PN=?;*TEMP2 
TEMPI  =  0. 

CALL  GAINCTFMP1  .FPL'.  TFMP2.VT.-1  .Fill 

PN  =  PN  ♦P.YFKK  **}!  /  CSORTCP.  *CPI  >*VT  )*TFMP2 
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PRC*PN*t<5PP*P'im 
PTC*P*J+f  '>TR  +  P7P> 

SNRAT=<PTT1  ♦PTC3/'<P\’R*PRC3 
361  SVDP=!P.*AL,’PinCSVPAT> 

365  TSM=F0StP1) 

P3D=0.‘ 

P3R=0* 

IFCSNDR.GF.TSM3P3D=l  . 

IFCP3D.LE.fi.  3  GO  TO  400 
370  CONTJVIIF 

XNR*AFIMt  CDTX,PTY3/DGX 
TFKP1  =  CX\'R-3.?3**J>/1  1  • 

IFCTEFP1 -LT. 1*0. 3  GO  TO  390 
P3R  =  t • 0 
•  GO  TO  400 

390  P3R=1 .-EXPC-TEFP1 3 
400  AT=0GX*CGY 

AT* AT  *  CNT+  CMT- 1 3+TSPACE/DGX) 

DSX=F0TS  <2  3 
DS Y=FDTS  <3  3 

XPAX=H/TAV<PHlD-TPFTAV/2. 3 
XM1N=H/TAN(PMID*TKFTAV/?. > 
YPAX=P.-*XMAX*TA«JCTHFTAW/?.  > 

YPIN=2  *  *  XM 1IM  +  TAM  <  THFT  AH/P •  3 
IFCCOSX.CT.D. 3.0R. CPSY.GT.O- 33  GO  TO  450 
AS  =  <XFAX-XMIN3*CCYr’AX-YFJ\’3/2.  3 
GO  TO  500 
450  DDX  =  XNA>;-xm|N 

IFCDSX.CT.O. 3  DDXsAFINl <DSX,DDX3 
DDY=<YFAX+YFTN3/4.  -Y 

IFCDSY.CT.O.  3  DDY=AVI\'l  C  CDSY/P. 3 ,PDY3 
IFC<DDY+Y3.G7. CYFIM/2. 33  GO  TO  A70 

as=?.*dpy*pdx 
GO  TO  500 

470  AA=<<DDY*Y3-CYMIN/P.  3  3 /TA\' CTHET AH/P.  3 
BB=AA*TA\'CTHFTAW/?.  ) 

AS=D0Y*D0X-.  5  »AA-*PO 
DDY=<YWX*YMl'J3/’4.  ♦  Y 

IF  CDS  Y.  (-T  .0.  5  DDY=AFIN 1 < CDSY/P.  3.DDY3 
IFCCDDY-Y3.GT.  CYMI.M/2.  33  GO  TO  490 
AS*AS  +  D!iY*DSX 
GO  TO  500 

490  AA=<  <PDY-Y3 -  C YFIN/P. 3 3/TAMCTHETAH/2. 3 
Bn  =  AA*TA\l<THFTAH/?. 3 
AS*A5+Df)Y*DDX-.  5*AA»PB 
500  ARAT10=AT/AS 

TEPP1  =  <700. /G)+ARATIO*DFLTAT 
IFtTFFPl.LT .180. 3  CO  TO  510 
P2*l .0 
GO  TO  5P0 
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510  P?= l.-FXP< -TEMPI > 

520  X*X/C2 
Y=Y/C2 
H=H/C2 
SR=SR/C2 
RFTURN 
END 

SUBROUTINE  INTFPP<TARGX,TAPGY,X,  Y. TAB, NX. NY.  AMS.  IfLaG.N) 

* 

PERFORMS  P 1 VAR  I ATE , INTERPOLATION 

DIMENSION  TAB  (N#  l)>X(l)iY(U 
I FLAG =0 

IFfTARGX-X<I )}  ?,3»3 

2  IFLAC=- 1 
AR0X*X< 1 ) 

CO  T0  A 

3  ARGX=TARGX 

4  lFtTARGY-YCtn  5,6.6 

5  IFLAG=-t 
ARGY=Y  < 1 1 
G0  TO  7 

6  argy=taroy 

7  DO  B  11=2. NX 
IF<ARCX-XX<mi  9»o,B 

B  CONTINUE 
I  I  =NX 
I FLAG = I 
ARGXaXtU) 

9  I  *  1 1 

DO  10  JJ=2.NY 
1F<ARGY-YCJJ) >  11,1!, 10 

10  CONTINUE 
JJ=NY 
1FLAC-=1 
ARGY=Y< JJ) 

11  J=JJ 

DY=Y(J>-Y< J-l 1 
CYI=-<ARGY-Y( J> 5/DY 
CY2=<ARCY-Y< J- 1 > )/DY 

dx=xc 1 7-xci  —  t  > 

CX 1  a  -  (ARC- X*  X  (  I  >  1/DX 
CX?  =  {ARGX-XCI-l Jl/OX 

ANS=CYI*CCX1 *TABC I  - t , J-l  ) -CX?*TAB C I , J- 1) > 

1  ♦CY?*<CX1*TAPC I-l, J)»CX2*TAP(I,J>> 

RETURN 

END 

SUBROUTINE  GA1N(X1 .XP.Y.VT, ITYPF.F1) 

OX=<X?-XI )/100. 

Y»0. 


V 1 1 -31 


x»xj 

00  40  I* 1,1 01 
IFCITYPF)  5,5»A 

5  TFMP=-X*X-..f>/C  VT*VT) 

T£YP=FXPCTFmP) 

IF C IT  YPF .LT.O)  TFMP=TFMP*X+*C?.*FI> 

00  TO  7. 

6  TFFP=SORT<C0SCX) )/CSINCX)**2> 

7  IFCI.FC.  !  )  00  TO  .70 
IFCI.FO.  ! 01  >  GO  TO  30 
*F<CI/t>*P-IJ  10,20.  30 

10  TEMP=2.*T!!>P 
00  TO  30 
20  TEMP=4.*TFMP 
30  Y*Y+TFMP 
40  X=X*DX 
Y=DX*Y/3. 

RETURN 

END 

SUBR0UT INF  SAPCFD0. IDT.FDT.FDB.FDFVT, 

*IDS,FDS, P2,P3D,P3R> 

C0f'M0N/T.ABLFS<'FOTAPC7PS) 

COMMON XTAOSFT/  ICTC50) 

C0MM0N/PL0C'<3/X>  Y.XY.P.SR,  5P»  AVPL0  »  ANPLT  >  OFLTAT 
DIMENSION  FP0C20), I 0TC20 ) , FDTC 20 > , FDPC20 ) , 
#FDFVT  C20) » !DSC20),FDSC25) 

CAS=FOS  C  1  9) 

ISFAS=IDSn> 

IBAND= 1 
FKT=J. 

CPi=?*A$iNn .  > 

C?*.30*R 
C3=I.6R9 
C8=. 01745 
C*3.F« 

DX=FDSC1  ) 

TAUPW  =  FRS  C  J  3  ) 

FLAMDA=FD$C9> 

IFCFLAMOA.LT. .025)  I BAND  =2 
IFCFLAMOA.LT. .0125)  IBAND  =3 
XXXX=0. 

DTX=FOT  C 1 >*C2 
DTY=FDTC2)*C2 
0ZFR  =  Ff)SC2) 

THFTAV=FDS  C 1 2 )*CB 
SLI =F0S  C20) 

SL2=FDS (21 ) 

SL3=F0SC?2> 

PX=FDSC5) 

TAUPI=FDSC14) 


VI 1-32 
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TAUP?=rr>s<i5) 

TSNsFDS (7  > 

PNR=FDS(3  ) 

PRF=FOS(A) 

PS I=FDT ( 1 1  )*CB 
RC=FDS(A) 

Z  =  63((7A50. 

AH=FDS(1R)*CB 

V01=FD3(S) 

H=FDQ (?)  *C2 
NT  =  IDT (  l  ) 

RN JC  =NT 

TSPACF=FDT  <1 0)*C2 
Y=FDO(l )+C? 

V.'S=FDS  (23)*(>076.*C2 
G=FDEVT(A) 

IF<Y.LT. Cl .5*WS>>  F0  TO  20 
C-0  TO  30 
20  Y=  1 .  S  +WS 
PRINT  25 

25  FORMAT  (  *  Y  TOO  SHALL*  INCRFASFD  T0  BE  J.5  SWATH  WIDTHS  *> 
30  0 AHHA=Y/Z 

222  SR=SCPT(?*Z*(Z+H>*( I -C3S f C-AHMA ) >*H+H> 

GRAr SP*T  AN ( AH ) 

SR=SPRT(SR*SR«-GPA*GRA> 

H0RIZ=SORTC (Z*H)**P-Z*Z) 

PHIH=ASIN(HORIZ/(Z+H>) 

PHI=ACOSC ( CZ+H>**2-SR**2-;>Z>/C-2*Z*5R>> 

PHI=PHI-CPI/2. 

FLP=C*TAUPW*.5/RC 
DC-Y=FLP/(C0S(PH1  )*C0S(AH>) 

D(?x=nx/cns(AH) 

DPMIN=AMINI  (DGX.OGY) 

RN JC 1 rDCHIN/CNT*(TSPACF*OTX>) 

RATS  I Or Af  IN  1  ( RNJC  *  PN,_  C  J  ) 

SICFAC  =  AN'AX1  (1  .  »  PATS  1  G  ) 

XTE=AMIN1 (DGX.DTX) 

XTY=AHIN! (DGY.DTY) 

ATF  =  XTF  *XTY 
AG=PGX  kDG Y 

ARF =AMA X 1 < (AG -ATF ) ,0.  ) 

PH I M=PH I -TPFTAV 

IFCPHIM.LE.PHIH)  PHIM=PHTH 

COHP=CPI /2. -PHIM 

ITCPHIM.LF.PHJH)  GO  TO  51 

C9MP=CPI -AS  IN (SIN (CP  1 X2. -3HIM)  +  CZ*H)/Z) 

51  GAMj-AM  =  CP!/2.  ♦PHIM-C0KP 
YMAX=Z*PHIH 

RAMAXYrACOS( (SR**P-(Z*H>+*2-Z*Z)/C-P*Z*<Z+H)> ) 
XY=Z*GAMAXY/C? 


VI 1-33 
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PHIM)  *PW!-<  Y^^HS  ■)/?. 

PAMMAXsAC^SC  crp/w  *9-?  +  7.*Z)/  (  -?«7!*Z>  ) 
x*Z*PA,vf'AX/C’ 

PAF=PZFR*(S  Mt  PM1  "I  )/S  IN  C PW I ) ) **P *SPRT 
* (COS (PH  I  )/C”S(PW!  v]  )  ) /COSt  AH) 

VP*FDO<5) 

AG1«AC 

ABFMO. 

IF (ARF •FO . 0.  )  GO  TO  POO 
105  flXX  =  AyAX! COGX.OGY) 

0N=APIN1 fOf X.Opy ' 

iFtoxx-LF.yp )  ro  to  ?oo 
IF(ON.LF.VP)  r-o  TO  110 
Ap=yp*toxx*n\’)/p. 

G0  TO  IPO 

110  AP  =  ON*CDXX*l-.'P)/p. 

120  AL'AHAXl C (AP-AP),0.  ) 

FM!  . 

AlrAP+Fl+AL 

IFtATF.GT. AP)  PO  T0  130 
ABE1 =ATF-AP 
G0  TO  MO 
130  ATFL=ATF-AP 

ABE1  *A1  -FMATFL-AP 
MO  AP 1  *F 1 + AP 
200  AT  =DTX*DTY 
PH1=PHI/C8 
N*15 

M6+IBAND 

BANO=IBANO 

S  IGTF  =  FOT  < I >*t.3/FLAMDA)**.3 
TFPP5  =  APS  t  PS I-AV) 

S  1  GTF  =  1  0  •  *ALOp1  n  c  s  I  PTr  ) »pw I / l  5 . -S/ AS. *TFWP5/C8 
IFtTFrPS.C-F.  •  7  *>5  /  )  SirTE  =  SICTFi-5/45.*TFHPS/CR 
*-5/A5  .  *  <  90  .  -Tc-vp^/CR) 

SlOTFclO.-n  <S1CTF/|P.  ) 

210  TF.VP3=AP/AT 

1 F <TF*'P3  .  CT  •  1  •  )  CO  TO  PPO 
SlPTF=TFf-P3*S  1  CTF 
220  S  IGTF  =  SIPFaC*SIPTF 
S 1  PR  1 =0. 

JJeJCTtlO) 

NX  =  JCT  Cl  1  ) 

N Y= ICT (12) 

IX-JJ 
I  Y=  1 X+NX 
1 7.  =  1  Y*N  Y 

CALL  INTER? (PH  I AND, FDTARC IX), FOTAPC I YJ.FDTAFt 17) , 
5  NX/NY.S irni , XXN>  ,NX> 

51CBl»10.**(SICnl/10.  ) 


VI 1-34 


PHI*PHJ*CP 

S  I GRF  *  A°F  1  *S  I GP 1 *SIN(PWI 1 
SIG£=AG1  *5IfR|i.SI\'(PMI> 

CALL  ATTF.\)C  IPAMO.CAS.W. tSFAS.PHl.TAUAR) 
TAUAR*  I  ./aO.*  +  <TAUAP/lO.  >> 

0*0. 

TEMP3=0. 

IFCSLl )  270.2P0.P70 
270  TFMP3=I  ./C10.**(SL1  /MO.  >  ) 

2R0  IFCSLP)  P90, 300,290 

290  TEMP3=TFMP3*1  .  /  (1  0.  *+  t  SL?/  I  0.  )  > 

300  1FCSL3J  3I0.3P0.310 

310  TFPP3=TFmP3  +  1  ./( IO.**<SL3/10. i> 

320  0=P.*TFXP3 

IFfO.FO.O. >  0  =  .  1 

PDO=TAUPI  *PX*GAF/M.  *CPI*SP*SR>*TAUAR 
PRTT=TAUAP+S IGTF*PDO 
•  PRPT  =T  AUAR-*S  T  0CF  *pr>0 
PRRPrTAHAR*S  IP(:*POO 

TFMPA=TAUPP*TAUAP*r-AF  +  FLAMDA**2/(A^CPI’*SR>** 

PATT=PRTT*TFMP4 

PART  =PPPT*TFMP4 

PARP*PPRP*TFMP4 

n_  =  .yt4*Fi  AMOA  +  SR/Dr-X 

FN=FL+PF  CC2*C3*VGl) 

PRT=FV*PAPT 

PTT=FW*PATT 

PRB=FN*PARB 

PSL=0*PRB 

SNRAT  =  C  PTT  +  PnT  )  /  C  PR9  +  PSL-*-PNR  ) 
SNRAT=IO.+AL0G1OCSMRAT> 

Y=Y/C2 

ANGL0=ATAMCY/X> 

SR=SR/C? 

H=H/C2 
SP=VG1*C3 
ANGLT  =  PH  I 
P3D=0. 

P3R=0. 

IFCSVRAT.GF.TSN>  P30= 1  •  0 
IFCP30.LF.0. )  GO  TO  /OO 
XNR=AMIN1 (DGX, Of- Y  >  / SORT  C  DGX*DGY > 

TEMPI =<XNR-3.2)**?/ 1  I . 

IF (TEMPI .LT. IPO. )  G0  T3  390 
P3R»1 .0 
GO  TO  A 00 

390  P3R=1.-FXP(-TFMPI  > 

400  AT=OGXtOGY 

AT=AT*CNT+(NT-I ) ♦TSPF CF/OGX > 

ASPECT  =FDS  <  1  7 ) 


WSM*VS 

AS  ■V'SH*VSH*  ASPECT 
500  ARATI0«AT/AS 

TV1FV«ASPFCT*  ■’Si*/  (SP»CP) 
TEVplc(700./n*ARAT!!1*TVlFV' 

IFITFHPt  .UT- I  ?0.  )  C-0  T0  510 

PP«l.O 

RFTURN 

510  P?  =  l #-FXPC-TEvPI  ) 

RFTURN 

END 

SUBROUTINE  ATTFNM  band.  CAS#  H,  I SFAS# PH I . ATT > 
C0MHON/T  AALFS''FDTA°(7?S) 

C0MHON/T ASSET' ICT (SO) 

•  JJ1 * ICT< 1 3 ) 

JJP* ICT  <  1  A) 

JJ3»ICTC I  9) 

JJA* ICT ( ?? ) 

JJS* ICT <  PS  > 

JJ6-ICTCPR) 

NXl * ICT <  t  4  ) 

NY1 * ICT  (IS) 

NX?=ICT< 17) 

NYP=ICTOR) 

NX3 - ICT  <  PO) 

N Y3* ICT  < ’ 1 ) 

NX4*ICI  <?3> 

NY4= ICT  CPA) 

NX5=1CTC?A) 

NYS* ICT  C  ?7 ) 

NX6=ICTC?9) 

NY6= ICT C  30 ) 

IZ=JJ?+NX°  *NY° 

1 1 =I2*NYP+ ISFAS 
IP*  12+ ISFAS 
CP Jrp. *AS  IMC  I  .  ) 

IF  CH.LT • FOTAP (11))  CO  TO  50 
1FCH.LT . FOTA” C IP) )  no  TO  70 

0FLS=Sl ( FOTAB  C II  )»  M. PH  I >-S  ! CPDTAPC IP),H»PHI ) 
00  TO  100 
50  DFLSrO. 

00  TO  100 

70  PFLS*$1 CFOTAPCI I ).H.PHI ) 

100  SUP=0# 


17.=JJ4*NXa«NY4 
00  POO  1*1 ,NY< 

N«=I2«.NY4*C T  SF  A  S -  1  >+ I 
SUH=SUV* FOTAP (N) 
IF(SUM.OF.CAS)  00  TO  P50 
POO  CONTINUF 
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P50  C«.OA 

IF  c  i pano . fo  .  p  i  r*.p 
IF  C  IPAVP. FP . 3  1  C.3.4 
U«JJ3»NY3*NX:i«NY3*C  1SFAR-'.  1*1 
ATT*C*FM  AP  C  1  I  1  ♦ . no  1  M1FLS*  f  I  -  I  >/R 

ifch.lt. fotaa <  ip)  i  ro  to  ?.-o 

DPLS  »S I CO.H.PHl 5 -SI CFDTAPC  I  I  l.H.PHI > 

C-0  TO  ?HO 

2A0  DFLS  *S 1  CO.H.PHl  1 
280  SUP=0. 

KC'O 

IFCOFLS.LT. 10000. >  KC*I 
IZ=JJS«NXS«NYS 
IFCNC.FO.il  I  ?  =  J JA  *NX6  +  N YA 
IFCPC.FO.il  NYS'NYA 
IFCPC.FO.il  NXSsNXf 
00  300  1*1 »NY5 
N*IZ*NYS«C ISFAS-t  1*  I 
SUK=SUP*FDTARCM) 

IFCSUP.CF.CAS1  CO  TO  350 
300  CONTIMCF 
350  IZ=JJ1 *NX| »MY1 

ILOC *  I  Z*.N  Yl  *(  1  HAND- I  1*  I 

ATR=FDTA«<  IL''C)*PFLS 

att>-att*atr 

RFTURN 

END 

function  si  chi .H.PHI ) 

CP  1 ♦  AS  INC  1 .  1 
Z«-6367A50. 

Si e-CZ«Hl*COSCCPI/:». -PHI l+SORTC C C Z  +  H>*COS CCP I /? . -PHI  1 >**? 
I* CZ«H1*»£>-CZ*HI  >+*2) 

RFTURN 

FNO 
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